Vortex Molecules in Spinor Condensates 
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Condensates of atoms with spins can have vortices of several types; these are related to the sym- 
metry group of the atoms' ground state. We discuss how, when a condensate is placed in a small 
magnetic field that breaks the spin symmetry, these vortices may form bound states. Using sym- 
metry classification of vortex-charge and rough estimates for vortex interactions, one can show that 
some configurations that are stable at zero temperature can decay at finite temperatures by crossing 
over energy barriers. Our focus is cyclic spin 2 condensates, which have tetrahedral symmetry. 
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In an image of a nematic liquid crystal by polarized 
light, one can identify topological defects of various topo- 
logical charges (see Ref. ^Ij). Bose condensates (see the 
books 0, 3) are starting to provide another context for 
studying topological defects: in the texture formed by 
the phase and spin of a spinor condensate [J, [a, [y]. (See 
also [7, 8] for reviews of the theory and experimental 
techniques being applied to spinor condensates.) 

A topological defect in the phase of a superfiuid is 
a quantized vortex. The discontinuity in the phase as 
the defect is encircled must be 27m for an integer n, 
and the circulation of the vortex is then n— . In a 

m 

single-component superfiuid, multiply quantized vortices 
(|n| > 1) are usually not stable. The widely known ex- 
planation is that the energy of a vortex is proportional to 
ri^. Thus a doubly quantized vortex {n? — 4) can lower 
its energy by splitting into two singly quantized vortices. 
Similar arguments can be formulated for multicompo- 
nent condensates, but we will find that some vortices in 
these condensates can be very long-lived in spite of hav- 
ing large energies. Our predictions are about condensates 
of atoms with spin in which the rotational symmetry has 
been weakly disrupted by a small parameter g, such as 
the interaction between the spins and the magnetic field. 
(We focus on the cyclic condensates of spin 2 atoms, see 
Ref. [9J.) 

Long-lived multiply quantized vortices are particular 
examples of composite vortices, which are made up of 
several vortices bound together. Ground states with 
complicated symmetries have many types of vortices [lOJ. 
Symmetry violating fields provide a force that can bind 
some groups types vortices together so that they form 
a "composite core" for a larger vortex. The compos- 
ite vortex might have two quanta of circulation (or have 
some other extra-large topological charge) . Such a vortex 
would be surprising if observed under a low resolution, 
because the smaller vortices that make it up would be 
hidden in its core; the behavior of the order parameter 
at infinity is determined entirely by the net topological 
charge. 

Some earlier theories about vortices in multicompo- 
nent condensates also describe vortices with asymmetri- 
cal cores which can be regarded as collections of closely 
spaced vortices. The parameter q makes this interpreta- 
tion even more meaningful: when q is very small, we find 



that the component vortices move far apart, so that their 
cores do not overlap, while still remaining bound. Ref. 
[ll| describes vortices that can occur in a condensate of 
two atomic states when there is an RF-field producing 
coherent transitions between the states. The binding of 
these vortices also comes from an asymmetry, but the 
asymmetry comes from the dependence of the interac- 
tion strength on the internal states of the atorns, rather 
than from an external magnetic field. Ref. [T^l stud- 
ies vortices of spinor atoms in a magnetic field, like us, 
but focuses on spin 1, and describes composite vortices 
as well. Since the scenario involves rotation as well as 
a magnetic field, the vortices would be held close to the 
axis of the condensate by rotational confinement without 
the magnetic field. With just a magnetic field, we show 
that the vortices are attracted to one another. 

Vortices and bound states are not hard to picture, by 
taking advantage of the fact that the state of a spinor 
atom can be represented by a geometrical figure. The 
appropriate shape depends on the type of condensate. 
For a ferromagnetic condensate, a stake pointing in the 
direction of the magnetization could represent the lo- 
cal state of the condensate. Other condensates can be 
represented by more complicated shapes. Now imagine 
a plane filled with identical shapes (tetrahedra, for the 
cyclic phase) , with orientations varying continuously as a 
function of position. This shape field (or "spin texture" ) 
together with a phase field would represent a nonuniform 
state of a condensate. If the shapes rotate around a fixed 
symmetry axis as some point is encircled then the spin 
texture has a topological defect at this point. Such con- 
figurations generalize vortices, because some of them are 
accompanied by persistent spin or charge currents. For 
each symmetry of the tetrahedron, there will be a vortex 
when the Hamiltonian is SU2 symmetric. (Each discrete 
subgroup of SU2 describes the vortices of some phase for 
atoms of some spin|13|; to find vortices for more even 
more complicated groups like SO^ or SOj, one might 
want to study gases of spin | atoms [14i,il5].) The vortex 
spectrum is decimated when a magnetic field is applied, 
since the field favors a particular orientation of the order 
parameter -e.g., tetrahedra may want to have an order- 
three axis aligned with the magnetic field. The spec- 
trum of vortices is then reduced from the full spectrum 
of "tetrahedral" vortices (based on arbitrary symmetries 



of the tetrahedron) to "field-aUgned" vortices, where the 
tetrahedra must rotate around this order three axis, so as 
not to lose their aUgnment with the magnetic field. We 
may introduce the ground state space Mq (of tetrahedra 
aligned with the magnetic field), where the interaction 
energy is V = Vmin, and the space of ground states of the 
SU2 invariant part of the Hamiltonian, M (of tetrahedra 
with arbitrary orientations) . A set of vortices can be as- 
signed a topological charge based on the loop traced out 
in one of these spaces by the values of the order param- 
eter on a circle around the set. Specifically, the topolog- 
ical charge describes the symmetry transformation that 
brings the order parameter back to its initial value as the 
set of vortices is circumambulated. (The relation between 
topology and symmetry is presented in [lO|.) From the 
general point of view, the reduction of the charge types 
from the tetrahedral to the field-aligned ones when q is 
turned on results because the smaller ground state space 
at nonzero q, Aiq, has fewer closed loops[52]. 

If q is very small, the excess energy of an order pa- 
rameter in A4 rather than Aiq is small. There is then 
a hierarchy of spaces Aiq, Ai and H (the whole Hilbert 
space, corresponding to arbitrarily distorted tetrahedra) 
with increasing energy scales. When g ^ 0, a vortex 
has to have a field-aligned charge at infinity because the 
tetrahedra must eventually move into A4q to avoid too 
big of an energy cost. But if q is small, the hierarchy of 
the order parameter space is reflected in the fact that a 
vortex with a field-aligned charge can have a composite 
core that may contain tetrahedral vortices, as in Fig. (TJ 
the core is sort of like the pulp and seeds of a fruit. Inside 
the fruit is a texture of arbitrarily oriented tetrahedra, 
almost as if q were equal to zero. The seeds are then 
the cores of tetrahedral vortices and the pulp is qualita- 
tively the same as the texture that would surround these 
vortices in the absence of the magnetic field. The net 
charge of the vortices has to be field-aligned so that the 
loop of order parameters that surrounds the whole core 
can move into Mq as r —^ 00. This picture becomes 
more accurate as q becomes small, since the vortex cores 
within the composite core are far apart in that case- the 
size of a region where the order parameter is in A^ — A^^ 
or 7i — A^ is inversely proportional to the energy scale 
for each space. As g ^ 0, the binding of the tetrahedral 
charges becomes weaker and weaker, until they become 
free from each other; each vortex therefore is described 
by independent degrees of freedom when q — 0. 

The optimal size Lq of a composite vortex results from 
competition between two forces-a confinement force from 
the anisotropy term and the familiar logarithmic inter- 
action of vortices. The symmetry breaking energy fa- 
vors minimizing the area over which the order parameter 
leaves Aiq, pushing the component vortices toward one 
another. On the other hand, it cannot compress them to 
a point since the Coulomb-like repulsion one expects of 
vortices keeps them apart. 

Some of these vortex molecules will turn out to be 
metastable. Ref. [121 mentions an interesting clue to 




FIG. 1: A composite vortex reflects the hierarchy of the or- 
der parameter space. In the white, grey, and black regions 
the order parameter moves from A4q to 7V( to TC which has 
the highest energy of all. The charges of the subvortices are 
tetrahedral charges, represented by F and the charges of the 
composite vortex is a "field-aligned" charge. 



such a phenomenon; namely there are multiple steady 
state wave functions describing a condensate with a given 
magnetization and rotational frequency; these local min- 
ima of the energy function can maybe be analyzed using 
the group theoretic metastability conditions we discuss 
in Section FlII CI For a spinor condensate, wave functions 
for states besides the ground state are experimentally 
important, since the experiments of Ref. [l6|, as well as 
the liquid crystal experiments of Ref. [17| , reveal compli- 
cated textures produced by chance; an initial fluctuation 
around a uniform excited state becomes unstable and 
evolves into an intricate nonequilihrium texture. So it is 
useful to analyze spin textures which are only local min- 
ima of the Gross-Pitaevskii energy functional (like the 
metastable vortex molecules considered here) as well as 
unstable equilibria (which take a long time to fluctuate 
out of their initial configuration). Examples of unstable 
equilibrium are described for single-component conden- 
sates by [l8|. The process by which textures form out 
of uniform initial states has been discussed in theoretical 
articles, including [19|, 1201 (on the statistics of the spin 
fluctuations and vortices that are produced from this ran- 
dom evolution) , [i21] (on the spectrum of instabilities) and 
[23 . |23| on the dynamics of spinor condensates. The ex- 
periments described in [2^ show that the patterns that 
evolve in rubidium condensates are probably affected by 
dipole-dipole interactions, though we are not considering 
these. Dipole-dipole interactions lead to antiferromag- 
nctic phases ;25|, which maybe can be described as ground 
state configurations of vortices. 

Besides just hoping for unusual types of vortices to 





FIG. 2: Geometrical representation of the tetrahedral super- 
fluid phases of spin 2 atoms. Two orientations of a tetra- 
hedron, corresponding to the spin roots of the ground state 
spinor, are illustrated and the symmetry axes are labelled. 
(The orientation of the coordinate a:xes, shown on the side, is 
the same in both cases.) a) An orientation with the vertical 
2-axis along an order three axis, corresponding to the ground 
state y/noX3 (see Eq. ((Ill), b) An orientation with the z-axis 
along an order two axis, corresponding to the ground state 
\/noX2 (see Eq. ([29}). For a magnetic field along the z axis, 
orientation a) is preferred when c > 4 and orientation b) is 
preferred when c < 4. Note that the three order 2 axes, la- 
belled A, B, C can be used as a set of three orthogonal body 
axes. 



form, one can make a vortex lattice by rotating a con- 
densate. The effects of rotating on a spinor condensate 
have been investigated theoretically in ^0| as well as the 
review [27| . Experiments can also make a single vortex of 
a prescribed type [2^, [2^, [SO]- Excited by these possibil- 
ities, physicists have come up with several types of vor- 
tices and topological defects they would be interested in 
seeing: skvrmions[3lL [32| . monopoles|33l|. textures whose 



order parameter- field lines make linked loops [34j, as well 
as the noncommutative vortices of the cyclic phase that 
we will be expanding on here[35|. 

The key to our discussion of vortices will be 
a geometrical representation of the order parame- 
ter, allowing us to visualize a texture of the cyclic 
phase as a field of tetrahedra with different orien- 
tations. (See Fig. [21) Without this represen- 
tation, a spin texture would be given by a spinor 
field {tp2{,x, y), '(/'i(a;, y),il^o{x, y), V-ila^, v), tp-2{x, y)Y] 
the fact that this spinor lies in the ground state manifold 
M would have to be described by a set of polynomial 
relations between the five components. A more reveal- 
ing way to represent a spinor is to draw a geometrical 
figure consisting of "spin-roots" (as in [Sa]), and in the 
cyclic phase, these spin roots form the vertices of a tetra- 
hedron. (A similar construction can be used to classify 
vortices in condensates of spin 3 atoms, see Ref. (37i].) 
Even without using the spin-root interpretation, one can 
justify using tetrahedra to represent order parameters in 
the cyclic phase because they are a concrete way of rep- 



resenting the symmetry of this phase. Ref. [35| worked 
out the symmetry group of a state in the cyclic phase by 
finding all the pairs n, a such that 
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/./I\ 



X3 



(2) 



V / 



The spinor y^noxs is a representative cyclic state, if no 
is the density of the condensate. These symmetries are 
the same as the symmetries of a tetrahedron oriented 
as in Fig. [2^; hence we may represent the state y^noXs 
by this tetrahedron. Any other ground state should be 
represented by a tetrahedron oriented so that its symme- 
try axes and symmetry axes of the spinor coincide. The 
appropriate orientation of the tetrahedron for a given 
ground state can be determined in an automatic way by 
calculating the spin roots. 

The Hamiltonian for spin 2 atoms in a magnetic 
field, H = JJ d^u[£^Vi;^ViP + VtotW] has a simple 
expression[9| in terms of the density n = V'^V'j nragneti- 
zation m = ■!/;^F'0, and singlet-pair amplitude 9 = tpl^P- 
(ipt stands for the time reversal of ^.) 

Vtotii') = i(an2-K/3m2 + c/3|0n-#ti^,V-M^V (3) 

The first three terms describe the rotationally symmet- 
ric interactions of pairs of atoms. The first one de- 
scribes repulsion between a pair of atoms and the next 
two terms describe additional, smaller interactions, that 
depend on the spin states of the two colliding atoms. 
These terms determine the spinor ground state in the 
absence of a magnetic field ^]. (The spin-dependent in- 
teraction strengths /?, c/? can be expressed in terms of 
the scattering lengths.) The properties of spin 2 atoms, 
which are described by this Hamiltonian, and of spin 
1 atoms have been investigated experimentally in Refs. 
|16l ISSl l39l . Ho, 14111: 4S | reviews more experimental phe- 
nomena. Ref. [33, H^ found values for a, /3 and c for 
^^Rb that are consistent with theoretical predictions, al- 
though even the sign of c is not known for sure, because 
c is small. The final term contains the chemical potential 

If f3 and c are positive, the ground state of a conden- 
sate of spin 2 atoms is cyclic. The deformation of a cyclic 
state due to a magnetic field is the simplest if we assume 
that a ^ |/?| and that c is on the order of 1. We there- 
fore assume c ~ 1, though c <C 1 for rubidium. The 
magnetic field influences the atoms through the fourth 
term of the energy, breaking the rotational symmetry of 
the Hamiltonian; this term describes the quadratic Zee- 
man shift due to a magnetic field B along the z-axis 



and q (X B^. (See Q for the explanation of why the 
quadratic Zeeman term q but not the hnear term is rel- 
evant if the condensate's initial magnetization is zero. A 
nonzero magnetization is described by a Lagrange multi- 
plier term —pJJ cPxip^ F^ip, which looks like a linear Zee- 
man coupling. We assume p — 0; a nonzero p should have 
similar consequences as a nonzero q, since both break the 
rotational symmetry.) 

Introducing a small q reduces the ground state space 
from M to Mq. This can be explained (see Sec. IIIAp 
by finding the modulation of the energy of an arbitrary 
tetrahedral state as a function of the orientation of the 
corresponding tetrahedron: 

^e// = (c - 4)-— 5 (cos'' ai -I- CDS'* ^2 + cos^ as), (4) 

where ai , a2 and a^ are the angles between the z-axis 
and three body-axes A, B, C fixed to the tetrahedron (see 
Fig. [2^) . The orientations that minimize this energy are 
the true ground states at nonzero q. If c > 4 (assumed 
until Sec. IIV|) . then the tetrahedron in its ground state is 
oriented with the z-axis perpendicular to a face as in Fig. 
[5^, with a ground state energy of Vmin = 4c0 ■ Thus 
the absolute ground state space Aiq = A4q3 contains 
all the wave functions that are arbitrary rotations about 
the z-axis (combined with rephasings) of ^/tiqXs- When 
c < 4, the ground states Aiq2 are rotations about z of the 
tetrahedron illustrated in Fig. ^p. In particular, when 
q y^ 0, there is a phase transition at c = 4, though there 
is nothing special about c = 4 in zero magnetic field. In 
this paper, we will mostly assume that c > 4 because our 
initial example occurs under this condition. 

Consider first a single vortex associated with the rota- 
tion through 180° of the ground state tetrahedron (see 

Fig. [2^) about the A axis, (w|,0, a/|)- Such a vortex 
is described at large distances by 



ip{r, (j>) 



^-4^iV2F. + F^] 



noX3 



(5) 



where r, (j) are polar coordinates centered on the core of 
the vortex. This vortex has an excess energy (relative 
to the ground state energy of the condensate) which di- 
verges with the condensate size. In fact, its Zeeman en- 
ergy density (with Vmin subtracted) is given according to 
Eq. O bylH 



(c-4)(g^) 
6/3c 



(sin^ 



(6) 



The tetrahedron is pointing in the wrong direction except 
along the positive and negative x-axis; and the integrated 
energy is proportional to the area: 



En 



q^R^ 



(7) 



where R is the condensate's radius. Such a vortex cannot 
be the only vortex in an infinite condensate! 
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FIG. 3: Illustration of the energy costs in a vortex molecule. 
The energy in the region occupied by the two vortices is dom- 
inated by the quadratic Zeeman cost and the energy outside 
it is dominated by kinetic energy. Increasing the molecule size 
decreases the kinetic energy and increases the Zeeman energy. 
The equilibrium size L is determined by minimizing the sum. 



The vortex described by Eq. ([5]) can form a partner- 
ship with another vortex of the same type, producing 
a molecule that can exist without costing too much en- 
ergy. This is because the net charge of the two vortices 
is compatible with the magnetic field. The combination 
of two 180° rotations of the tetrahedral order parameter 
about the A axis is a 360° rotation about the ^-axis, but 
since any rotation axis is a 360° symmetry of the tetrahe- 
dron, the rotation axis on circles of radius r enclosing the 
vortices can be tilted continuously relative to the tetra- 
hedron as r increases until it becomes the R axis instead. 
Then the tetrahedra align with the magnetic field far 
from the vortices and stay in Mqs- Each vortex screens 
the part of the other vortex's charge that produces the 
large Zeeman cost, as illustrated in Fig. [3l The Zeeman 
energy in the region around the vortices, where the tetra- 
hedra are still tilted, may be estimated by replacing the 
total condensate size R in Eq. ([7]) by the diameter L of 
this region. The Zeeman energy tries to pull the vortices 
toward one another but the elastic energy cost of rapid 
changes in the order parameter (the gradient term in the 
Hamiltonian) opposes this tendency: The tetrahedra ro- 
tate twice as fast around circles beyond L, where the two 
vortices act in concert. Therefore the elastic energy in- 
creases as L becomes smaller; this leads to the Coulomb 
repulsion between the vortices, 27r2i^^ In .^ — tt^^^i^ In — , 

^ ^ ni a ni ac ^ 

where m is the mass of the atoms in the condensate. The 
energy of the vortex molecule is therefore 



E^k 



q^ 

P 



nofr L 
n In — 

m Ur 



cnst. 



(8) 



where fc is a numerical constant. The equilibrium size 
can be determined by minimizing over L. The quadratic 
Zeeman force binds this molecule together while the re- 
pulsion keeps the vortices from merging. 

If the vortices were to coalesce, then they could react 
to form a set of vortices that are not bound by the Zee- 
man energy; one possible set of decay products is three 
vortices each involving a rotation through 120° about 



the R axis (which have the same net 360° rotation as the 
original pair of vortices). The vortex molecule of the two 
A rotations is metastable because thermal fluctuations 
may overcome their Coulomb repulsion and push them 
together, leading to such a fission process. 

The rest of this paper elaborates: it gives a qualita- 
tively correct expression for the spin texture surround- 
ing the molecule just described and determines how this 
molecule actually decays. In order to do this, we will 
give some more general results: first Section HI summa- 
rizes the non-commutative group theory of combining 
vortex charges and a classification of the tetrahedral vor- 
tices (the vortices that occur by themselves when q = 
and in clusters when q ^ 0); then Section HIl estimates 
the elastic energies and Zeeman energies of such clusters 
as functions of these vortex charges; finally Section IIIII 
gives criteria determining which types of tetrahedral vor- 
tices form bound states or metastable states (see Section 
IIIIC|) . The last two sections illustrate the criteria with a 
few additional surprising examples (see Section HV]) and 
give some basic ideas about how to observe metastable 
vortices (Section |V|. 



I. TOPOLOGICAL CHARGES 

Vortices are simplest to understand for an interaction 
energy that has a single manifold. A/", of ground states 
and no hierarchy. The order parameter must move into 
A/" far away from any spin texture. Vortices are classified 
by the topology of the circuit traced out in M by the 
order parameter on a large circle containing a vortex or 
set of vortices. Although there are many ways to add 
wiggles to a given circuit, only the topological structure 
is important, leading to a discrete set of possible vortex 
charges. Fig. [4] shows how a circuit may be tangled with 
holes in the space A/". As a vortex evolves, such a circuit 
can evolve only into other circuits that are tangled in the 
same way; thus the vortex charge is conserved. 

This generalizes circulation-conservation in a single- 
component condensate: The circulation quantum num- 
ber n for a set of vortices is also the number of times that 
the wave function at infinity winds around the circle that 
minimizes the Mexican hat potential|43|. The tangling 
in a multicomponent condensate can be described by a 
group of generalized winding numbers around the ground 
state manifold A/", tti (A/") , the "fundamental group" (see 
Fig. El). 

Besides the conservation of topological charge, two fur- 
ther properties of vortices follow from the geometry of the 
internal ground-state space. First of all, the net winding 
or topological "charge" of a set of vortices can be found 
by multiplying their charges together, using the definition 
of multiplication for the fundamental group. (In a gen- 
eral space, the fundamental group has a multiplication 
operation defined by splicing circuits together.) This rule 
is the generalization of adding the n's of the individual 
vortices in the scalar order-parameter case. (The funda- 




FIG. 4: Tlie order parameter at infinity of a vortex winds 
around a loop F in the order parameter space M. Realistic 
order parameter spaces are usually more symmetrical. 



mental group is noncommutative so one has to multiply 
the charges together in the right order, see Appendix [B]^ 
Secondly, the energy of a vortex can be estimated as a 
function of the winding behavior at infinity. For a vortex 
that minimizes this energy, the order parameter travels 
along a geodesic in A/". Since the interaction energy V 
is constant (and equal to its minimum Vmin) at infinity, 
the energy of a vortex is determined by the elastic en- 
ergy, that is, the cost of variations in t/i as a function 
of the azimuthal angle (j>. The closed loop traced out 
by ip{R cos (f), Rsmcf)) will relax to make this energy-cost 
small. When it shrinks as much as is possible without 
leaving A/", it becomes a geodesic in A/". The energy of 
the vortex is related to the length I of this geodesic by 



i^.(^ln^)i^ 
m ar 47r 



(9) 



where R and Oc are the radii of the condensate and the 
vortex core. (Note that replacing / by 27m gives the stan- 
dard expression for a vortex in a scalar condensate.) 

A small magnetic field introduces hierarchies into the 
order parameter space, leading to spin textures that wind 
around one manifold at an intermediate length scale and 
around a smaller manifold far away. 



A. Spin Textures and their Vortices 

Let us consider how nested vortices hold themselves 
together in a continuous "spin texture." We will try to 
give a general argument showing how topology and the 
energy-hierarchy of the order parameter spaces implies 
that any texture is made up of a set of composite vortices 
which in turn are made up of clusters of nearly point-like 
vortices. 

The order parameter subspaces in order of increas- 
ing energies are the two-dimensional AAq, the four- 
dimensional A^ and the AF -\- 2 dimensional H (for 



a) 




= 



= 27r 



b) 




FIG. 5: Classifying vortex topologies, a) The combined topol- 
ogy of a set of vortices can be found be seeing how the order 
parameter changes around a closed loop parameterized by <j)- 
The loop is drawn with a small gap between (jj — and (jj — 2-k 
to indicate that (?(0) and 0{(j)) vary continuously only between 
and 2tt. The topological charge is defined as the jump in 
the values of g and 9 across the gap. b) Vortices that cannot 
exist by themselves in the presence of a magnetic field can still 
sometimes appear in clusters. The charges Fi, F2,F3 in one of 
the clusters can involve rotations around arbitrary axes, but 
the combined charge Q has to use the direction of the mag- 
netic field for its rotational axis to avoid a large energy-cost. 



spin F), with corresponding energy scales V — Vmin — 
0, e{q) = ^ (see Eq. (g])) and /3ng (see Eq. (O; as long 
as the density of the condensate does not vary, /3 rather 
than a sets the energy scale). The condensate can move 
out of the ground state into one of the higher-energy 
subspaces if forced to by topology, spending less space 
in the manifolds with the greater energies. (This is like 
saving eggnog-with all that fat!-for one week out of the 
year.)For a particular subspace Mi in this sequence, in 
which V{ip) — Vmin ~ Ei, the size Li of the regions where 



h np 



IS 



-0 is in Ali is typically Li ^ \/j': where K 
the elasticity of the condensate. This relation can be 
guessed at without understanding anything about the 
field configuration in such a region; just assume equipar- 
tition between kinetic and interaction energies, so that 
ei is equal to ■^, a typical scale for the kinetic energy 

density, ^^iVf/'p. Consequently, the composite vortex 
cores in Fig. [I] where -0 is in A^ — A^^, have diameters 



(10) 



where the energy scale is taken from Eq. (|4]) and the 
component cores, where ip varies through Ti — A4, have 
typical diameters a^ ^ </-^. 

V "oP 

The more complete version of the equipartition ar- 
gument, given in Sec. IIIIBl estimates the total en- 
ergy of the vortex as follows: the kinetic energy out- 
side the composite core is about K-^luj^, according 
to Eq. ([9]), the magnetic field energy inside the core is 
about e{q)L^. Minimizing the sum, which has the form 

—ciK\iiL + e{q)L'^+C2, gives L ^ \/'Ma- The analogous 




argument also applies to vortices in a single-component 

condensate, giving the core size . /-^. Because the com- 

posite cores in a cyclic condensate have some substruc- 
ture (the component vortices), the actual coefficient ci 
of the logarithm for the cyclic condensate will turn out 

,2 

to be less than ^ once the kinetic energy of the vortices 
inside the core is included (see Section [TrBJ) . 

If q is very small, the component cores are much 
smaller than the composite vortices, and may be regarded 
as points. The wave function can therefore be approxi- 
mated by a field that always stays in A4 except at "sin- 
gularities" corresponding to the component cores: 



i^{x,y) 



^-ia{x.,y)n(x,y)-F ^ie{x,y) 

except at "points". 



"0X3 



(11) 



where xs is the spinor defined in Eq. ([2]). Over distances 
much greater than Lg even the composite cores seem to 
dwindle to points, suggesting the following approxima- 
tion: 

V'(a;,y)«e-^"(^'^)^^e*«(^^^)y^X3 

far from clusters. (12) 

These expressions generalize [4J| the "phase-only approx- 
imation" for ordinary superfluids[43|; for example, Eq. 
(fTTj) parameterizes elements of the space Ai in terms 
of the symmetries of the q = Hamiltonian, rotations 



ij 



-iaF-Ti 



ip, and rephasings, ip — > e ip, which are ap- 



plied to a representative state. Vortices are points which 
look like singularities at the level of resolution exposed by 
one of these approximations, although what one identifies 
as vortex singularity depends on which level of resolution 
one uses! 

Composite vortices result from the order parameter 
being forced out of the minimum-energy space Mq by 
topology: If the order parameter goes around a hole in 
A4q on some circle in the condensate, then inside the cir- 



cle, the order parameter has to leave A4q. If Eq. (fT2|) 
continued to hold all the way to the center of the circle, 
then there would be a singularity since the angular vari- 
ables a and 9 would run rapidly through multiples of 27r. 
So we can think of the circle as the core of an extended 
vortex in the field in Eq. (fT2|) . The singularity is actu- 
ally filled in by a field of tetrahedra with higher-energy 
orientations, described by Eq. (jlip , just as the wave func- 
tion of a vortex in a single-component condensate avoids 
singularities by vanishing in the core. In the fiUed-in 
region of freely-oriented tetrahedra, there are also some 
elementary "point" -vortices, actually spread out over the 
distance Uc- These result when the order parameter gets 
tangled around holes in M , and then has to move outside 
of this space. These cores are smaller than the compos- 
ite core because of the large energy scale associated with 
M . The fields of these vortices can be described qualita- 
tively by setting q = since kinetic energy surpasses the 
anisotropy energy very close to the cores. 

The cores of the component vortices are filled to al- 
most the same density no as the rest of the condensate 



(the spin-independent repulsion a favors uniformity since 
a ^ f3). Vortices with this property are often referred 
to as "coreless," but we use the word "core" in a more 
general way, so that every vortex can have one. Instead 
of defining core as a region where the density vanishes, 
our definition identifies a core as a region where the wave 
function departs from any particular form. (In the core 
of a "coreless" vortex in a ferromagnetic condensate, the 
order parameter becomes polar rather than ferromag- 
netic.) The core of the component tetrahedral vortices 
is the region where the order parameter leaves M. Us- 
ing the absolute ground state space Mq leads to a dif- 
ferent definition of vortex cores of field-aligned vortices: 
the core is where Eq. (fT2|) breaks down. This is why 
we call the region surrounding the tetrahedral vortices a 
"composite core." The ordinary definition of a core arises 
when we focus on another space, S, the sphere defined 
by |'0I — V^Q- (This space minimizes the largest of the 
interaction terms, ^an^ — nn.) From the perspective of 
S, a core would be a region where the wave function van- 
ishes to avoid being discontinuous; however, in a spinor 
condensate with a 3> /3, the wave function never has to 
vanish because a closed loop cannot get snagged in the 
surface of the simply-connected sphere S. 

Now we have to classify both the topologies of the 
point vortices (or "tetrahedral charges") and the com- 
posite vortices (or "aligned charges" ) , using the method 
described e.g. in [lO|. For the point vortices we must pa- 
rameterize the symmetry group that generates the com- 
posite core spin textures by a simply connected group G* . 
We will take G* = {(.g,6')|.g e SU2 and 9 e R} where 
multiplication is defined by (51, 0i)(52, ^2) — (<?i52,^i + 
02)- We parameterize G* as follows, 



Die-^ 



^lu — ian-. 



(13) 



allowing us to regard G* as a (redundant) description of 
the symmetries of the Hamiltonian. Using this represen- 
tation, one can define the net vortex charge inside any 
closed loop A (see Fig. [5^); one simply expresses the wave 
function along the loop in the form 



V(0)=i?(.9(0),e(0))V^X3 



(14) 



where (f> parameterizes the loop (0 < (p < 27r). In order 
for Eq. (fT4| to be continuous when the circuit is closed, 

i/;(0 = O) ^t{j{(j)^2Tr), or 

D{g{0),e{0))x3 = D{g{2n), 9i2n))x3. (15) 

It follows that the "classifying group element" or net 
"topological charge" inside the curve 



r(A) = (e 



— ("c-iaA- 



,0x) = {g{O)-'g{27T),9{27T)-e{Q)) 

(16) 
is a symmetry describing the net rotation and rephas- 
ing around the closed loop. (Let us also use the 

briefer notation g^ = e~*"^~2— for the net rota- 
tion.) Since the tetrahedron has only 24 symmetries 



(see Sec. II B|) . this makes for a tractable classifica- 
tion of the topologies. These topologies form the group 
of tetrahedral charges, multiplied together according to 
{gi,0i){g2,02)^{9i92,Oi + 02). 

When q is very large then Eq. P^ has to apply ba- 
sically everywhere and when q is small then Eq. (jl2|) 
has to apply outside of bound clusters of vortices. In 
order to classify the topology of vortices in the former 
case or of vortex clusters in the latter, we need a sim- 
ply connected group that parameterizes the q ^ sym- 
metry group e*^e-*"^-; we take G* = {{a,e)\a,e e R} 
and use the mapping D{a, 0) = e*^~*"^-. Then we write 
-0(0) = Dq[a{(j)),0{(i)))y/n^xs and define 



Q{\) = {ax,exh - (a(27r) - a(0), 0{2^) - 0(O))3 (17) 

where the subscript 3 indicates that the magnetic field fa- 
vors the 3-fold symmetric orientation illustrated in Fig. 
[2^ (since we are assuming c > 4), and the rotation 
through a\ is understood to be around the field axis. 
The possible values for Q can be referred to as "field- 
aligned charges" since they describe the topologies of 
vortex-fields in which the tetrahedra keep the orientation 
favored by the magnetic field. 

Now let us consider the form of the fields near a tetra- 
hedral vortex core. Near vortex z, the spin texture will 
be rotationally symmetric, and given by 



i, 



"0X0 



(18) 



where <f> is now the azimuthal angle cf) centered at this 
vortex. Because the tetrahedra near this vortex may be 
tilted, we use y^noXOi a generic member of the cyclic 
order parameter space. The vortex is azimuthally sym- 
metric since the parameters 0i,ai, h'^ are constants. The 
rotation axis n'^ is a local symmetry axis for the possibly 
tilted tetrahedra. To reduce the possible symmetries to a 
finite set, we should relate the spinor xo to the spinor xs 
corresponding to the tetrahedron as oriented in Fig. [5^. 
If we write xo = D{R,^)x3 for an appropriate rotation 
R and phase £,, the vortex in Eq. (flSl) can be written 



V'(0)=£'(i?,e)e'^'^e 



-i-^4>rii-Y 



X3, 



where 



fii = R-\n[) 



(19) 



(20) 



and we have used the transformation rule for angular 
momentum: 

3 
D{R,O^F^D{R,i)^Y.^^i^o- (21) 

(The right-hand side uses the SO3 matrix Rij associ- 
ated with the rotation R.) Continuity implies that the 
rotation axis n^ is one of the finitely many symmetry 
axes illustrated in Fig. [2^; furthermore, according to the 
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above scheme, the group element that classifies this vor- 
tex is Ti — (e~*"'^2^, e*^'). There are only a discrete 
set of possible charges when we use fii, the axis relative 
to the body axes of the tetrahedron rather than n^, the 
axis relative to the lab coordinates. 

Eq. P^ expresses the vortex as a product of a con- 
stant matrix D{R,^) (the phase ^ is unimportant) and 
a standardized vortex configuration. The transformation 
D(R,^) rotates the standardized configuration in spin 
space, changing both the rotation axis (from hi to the 
local axis n^) and the orientation of the tetrahedra. 

In the distant surroundings of a cluster of vortices, 
the spin texture will again have a rotationally symmetric 
form. The group element describing the change in the 
order parameter as one tours the loop A enclosing the 
entire cluster is given by 



r(A) = nr. 



(22) 



This algebraic law has a few consequences. First, it leads 
to a conservation law that constrains vortex alchemy: the 
net charge r(A) has to be conserved as the vortices inside 
the loop combine and metamorphose, since the topology 
on the loop cannot change suddenly. (There is only a 
discrete set of vortex charges because charges are defined 
using the body coordinates.) Second, Eq. ([^ restricts 
the types of vortices which can form a cluster when q ^ 0. 
(See Fig. [51d.) Outside the cluster, the tetrahedra must 
be aligned with the field (see Eq. ([H]), so the charge 
is described by a field-aligned order parameter Q{X) — 
{a, 9)3. Eq. dig) requires that 



n^^ 



(23) 



and J2i ^i = ^(mod 27r). The group elements must multi- 
ply to form a rotation about the z-axis to avoid the large 
Zeeman cost outside the clusters, or "composite cores" 
of the vortex molecules. 

Let us now review the example in the previous sec- 
tion: the vortex molecule was made out of two vor- 

■ ' -AV2a^+a,) r,\. „:„„„ -n2 _ 



(e" 



, 0); since F^ 



tices of type Fi 

(e"" '•3^"%0) = {-id,0) = (e-"'"%0), the spin tex- 
ture can align with the magnetic field outside the pair 
of vortices and Q = {2tt, 0)3. This uses the fact that, in 
SU2, all 27r rotations are equal to —id, where id is the 
2x2 identity matrix. 



B. Notation for Vortices wth and wthout a 
Magnetic Field 

Let us first assign names to the tetrahedral charges 
corresponding to closed loops in A^. At g = or within 
vortex clusters where q can be neglected, the topological 
charges are described by a pair F = {g,d) (defined in Eq. 
dHl)) where g = e''"^ € SU2 and 6* e M. Note that 



for the phase of the order parameter to be continuous, as 
Eq. P^ requires, the allowed values of 6 and g must be 
correlated [33, [S^, so only a discrete sequence of phases 
may accompany a given rotational symmetry. Also note 
that in SU2 , rotation angles are defined modulo 47r rather 
than 27r. The twelve symmetries of the tetrahedron ac- 
cording to the ordinary method of counting become 24 
because, e.g., a clockwise 240° rotation around an axis 
is distinguished from the counterclockwise 120° rotation. 
This is more than a technical point: the vortex where the 
order parameter rotates through —240° about an order 
three axis cannot deform continuously into one where the 
order parameter rotates through 120°, and it has more 
energy as well. On the other hand, the a = An "vortex" 
can relax continuously to a state free of vortices. The ne- 
cessity of using SU2 instead of SO3 is the biggest surprise 
to come out of the topological theory. 

Since g must be a symmetry of the tetrahedron cor- 
responding to ^yriQXa, illustrated in Fig. [^tt, we can 
describe g by indicating its rotation axis using the la- 
bels from the figure. We refer to the minimal rotation 
around a given axis using just the label of the axis; 
hence S, P,Q,R refer to the rotations through 120° coun- 
terclockwise as viewed from the tips of the correspond- 
ing arrows, and A, B, C refer to counterclockwise rota- 
tions through 180°, about A,B,C. Rotations through 
larger angles can be written as powers of these rotations. 
Therefore P^ is a 240° rotation; also P^ = A^ = -id 
since 360° rotations around any axis correspond to —id 
in SU2- We find it convenient to describe each rotation 
as a rotation through an angle a around some axis, where 
— 27r < a < 27r. (Positive and negative a's correspond 
to counterclockwise and clockwise rotations respectively.) 
An arbitrary rotation angle can be replaced by an angle in 
this interval using the fact that 720° rotations in SU{2) 
are equivalent to the identity. For example, the 480° 
counterclockwise rotation P^ is the same as the —240° 
clockwise rotation P^^ because P^, a rotation through 
two full turns, corresponds to the identity of SU{2). On 
the other hand, the 240° counterclockwise rotation P^ is 
not equivalent to the clockwise 120° rotation around the 
same axis, since the corresponding SU{2) matrices differ 
by a minus sign. 

Now we can list all the pairs of rotations and phases 
which are allowed by the continuity condit ion. Eg. ( j 1 5 p . 
The possible vortices according to Ref. [33, [33] are 
(i?™, ^^ -f 27rn) and (A, 27rm) where n and m are in- 
tegers, as well as the corresponding vortices with R re- 
placed by P, Q, or S* and A replaced by B or C. 

One can work out explicit expressions for the SU{2) 
elements corresponding to given rotations. For example, 
let us find the SU{2) element corresponding to A] since 
the rotation angle is 180°, 



A 



2 = la ■ cr. 



(24) 



where a is the ^-axis. Note that the tetrahedron 
has its vertices at (0,0, -1),(-^, 0, §), (^,^1, |), 



(-^, — i/l, i). The A axis bisects the segment connect- 
ing the last pair of points; the midpoint of these two 
points is 



2^^ 



^4»4) 



(25) 



The unit vector a is obtained by normahzing this vector, 
so 



a-(V|oV^)- 



Hence 



A 



-M/-a 



Vr- 



(26) 



(27) 



The net charge, Eq. ([22|) . of a set of vortices resuhs 
from multiplying the matrices g for the vortices of the 
set. The result can be identified as one of the rota- 
tions A",i?",P", etc. This procedure completely deter- 
mines the SU(2) product element, whereas the geometric 
method of applying the appropriate sequence of rotations 
to a tetrahedron does not determine the sign of the SU{2) 
matrix. 

We will describe a vortex reaction with the following 
notation, 



(_,rf,0)^(i?,y)*(i?,y)*(i?. 



4n 



(28) 



Each factor describes the charge of an individual vortex- 
atom, rather than a cluster of atoms, although {g, 9) can 
be used to describe the net charge of a set of vortices as 
well, as in Fig. [5^. This reaction describes onne point- 
vortex breaks up into three vortices. The * is just a sepa- 
rator between the different reaction products, reminding 
one to check conservation of charge by multiplying both 
sides of the reaction out. 

Another useful cyclic spinor is 



/ ^ \ 




X2 



V2 


1 

2 



(29) 



V ^ / 



This spinor corresponds to the tetrahedron in Fig. 
^jp. Its vertices are at the points of the form 
(±4=, ±4=, ±4=) if we restrict the choices of signs so 
that there are always or 2 minus signs. The fact that 
the A, B, and C axes of this tetrahedron correspond 
to the X, y and z coordinate vectors makes the spinor 
\/naX2 especially convenient for determining the conse- 
quences of the quadratic Zeeman term in the next sec- 
tion. The orientation of the tetrahedron in Fig. [2h 
is also very convenient for working out the group of 
charges, since the expressions for the symmetry axes 



are so simple. (E.g., p = (-i=, --i=, --i=), since ver- 
tex P is in the x>0,?;<0, z<0 octant. Thus 

p _ -if pa- _ l-i<7x+i<yy+i'^z \ 

At nonzero magnetic field, only rotations around the z- 
axis are symmetries. When c > 4, the ground state space 
is A4q3, consisting of rotations and rephasings of ^ynoX3^ 
as in Eq. p2)) . Vortices are described as in Eq. (fT7|) 
by an ordered pair {a, 6)3 describing the rotation and 
rephasing angles of the vortex. The subscript 3 is used 
to indicate that the z-axis is an order three symmetry of 
the c > 4 ground state tetrahedron. The continuity of 
the phase limits the vortex types to the form 

yTTTn TTi 

Q = («,^)3 = (^-,2^(y+n))3. (30) 

When c < 4, minimizing Eq. ^ implies that the 
magnetic field axis is an order two symmetry, and the 
ground state space is M.q2, the rotations and rephasings 
of ^/nox^i■■ Now {a, 6)2 specifies the vortex types. The 
possibilities are 



Q = (a, 9)2 = (7rm, 27111)2. 



II. ENERGIES AND SYMMETRIES 



(31) 



This section considers the effect of the magnetic field, 
which binds vortices, and the kinetic energy, which keeps 
the bound vortices from merging altogether. These are 
both included in the full energy function 



n 



d2j. v^tv^ + 14„,(^), 
Am 



(32) 



where Vtot is given by Eq. ^. When S 7^ 0, the Zee- 
man effect introduces an extra phase boundary dividing 
the cyclic phase into two phases, with a phase transition 
at c — 4 where the tetrahedron changes its orientation 
relative to the magnetic field. (This result applies for 
moderate magnetic fields; very low magnetic fields cause 
different transitions [26J.) 



The Anisotropy potential from the Quadratic 
Zeeman effect 



Let us determine which orientations of the tetrahedron 
will be preferred by a magnetic field along the z axis. The 
preferred orientation can be calculated for small q from 
an effective potential which is a function of the orienta- 
tion of the tetrahedron. As long as 



q < noP, 



(33) 



the tetrahedron will be only slightly deformed. It will 
move into a space M.' displaced by a distance on the 
order of — ^ from the space AA of arbitrarily oriented 
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perfect tetrahedra. The spinors in the distorted space 
are given by 



iIj = ^D{R,0x2+S^, 



(34) 



where D{R, ^) is the spin two rotation matrix correspond- 
ing to the rotation R of space, multiphed by a phase. The 
distortion 6ip depends on the orientation R. Eq. (fTTjl . 
which omits the deformation, is a harmless shorthand 
description, emphasizing the orientation of the tetrahe- 
dron. In this section, we use y^noX2 as the standard 
spinor orientation instead of y^noxs to simphfy calculat- 
ing the energy; conveniently, the body axes A^ B,C' of 
the corresponding tetrahedron for the former state are 
aligned with the coordinate axes x, y, z. The body axes 
of the rotated state D{R)y/noX2 (which make the an- 
gles tti, a2, as with the z-axis) are thus i?(x),i?(y), i?(z). 
Therefore the z component of the spin, in terms of the 
components of the spin along the body axes, is 

D{RyF,D{R) = cosaiF^ -fcosaaFj, -l-cosaaF;,. (35) 

At first order the quadratic Zeeman effect does not 
have any dependence on the orientation of the tetrahe- 
dron because a tetrahedral spinor is "pseudoisotropic," 



i.e., xlFiFjX2 



25, 



The first order energy is thus 



given by 

< qF^ > « noqxlD{R)^F^D{R)x2 

3 

~ noQ ^ cosaiC0sajxlFiFjX2 

~ 2no9, (36) 

which does not prefer any orientation of the tetrahedron. 
In the last step we used 



y cos ai 



(37) 



which follows from the fact that z, a unit vector, has 
body-coordinates (cosai, cos 012, cos as). 

To find the second order energy due to the quadratic 
Zccman effect, which will break the tie, we have to find 
the deformed state and its energy. The deformation 6'>p is 
determined by minimizing the total interaction and Zee- 
man energy in Eq. ^ for each given orientation R. Now 
if the deformation is not restricted somehow, the "defor- 
mation" which minimizes the energy will be very large, 
involving the tetrahedron rotating all the way to the ab- 
solute ground state. We therefore allow deformations only 
of the form 

6ij^dD{R,Ox2 
+ ^{aD{R, i)F^X2 + bD{R, 0FyX2 + cD{R, 0F^X2) 
+ {e + if)DiR,0x2u (38) 



where a, 6, c, d, e, / are real numbers. These terms cor- 
respond to the excitation modes found by [43]. This 
correction only perturbs "0 in 6 of the 10 directions in 
the Hilbert space. The other 4 directions are accounted 
for by the rotation R and the phase ^ which would be 
Goldstone modes when q ~ 0. (The energy remains ^- 
independent even when q ^ 0.) The particular six stiff 
deformations in Eq. (|38p are chosen because they are or- 
thogonal to infinitesimal rotations and rephasings of the 
tetrahedral state. We have to find the deformations that 
minimize Vtot-, Eq- ([3]), for each rotation R. 

To evaluate Vtot, note that D{R, ^) cancels from all the 
terms in the energy except for the Zeeman term, where 
one can use Eq. (l35|) . The resulting expression for the 
energy density reads 

3 

— q 2_\ cosai COS ajip^FiFjip (39) 

where V^ = y/nEX2 + '^X2 + ^ {aF^X2 + bFyX2 + cF^X2)] + 
(e -I- if)x2t is the perturbed wave function without the 
rotation. The effective potential Eq. ^ is obtained 
by minimizing Vtot over a, 5, . . . while keeping R fixed. 
Further details are in Appendix [A) (Working with ip 
is equivalent to fixing the orientation of the tetrahedron 
and rotating the magnetic field.) 

The effective potential suggests an analogue of the 
magnetic and charge healing length in condensates with- 
out magnetic fields (the "tetrahedron tipping length" )- 
when the tetrahedra are rotated out of the appropriate 
ground state orientation at the edge of a condensate, the 
energy in Eq. (J4| returns the order parameter to A4q 

within the distance L„ -- ^.F^. 

y q \ m 

Now the ground states can be found as a function of 
c. When c > 4, a short calculation shows that Eq. (JH) 
has its minimum at cosa^ = ±-7^, i = 1,2,3; i.e., when 
the magnetic field is along the line connecting a vertex 
of the tetrahedron to the opposite face or vertex. Hence 
the order parameter space Mq is as given in Eq. (|12p. 
When c < 4, the effective potential is minimized by an 
orientation in which the field points parallel to the line 
joining a pair of opposite edges (see Fig. [^[54|. 

The dependence of the ground-state orientation on 
c can be understood intuitively using the geometrical 
representation of Ref. [3(5]. In a tetrahedral state, 
m = ip'^'F'ij} = and 9 = tpjip = Q, minimizing the in- 
teraction energies in Eq. ([3]). When a magnetic field is 
applied, the base of the tetrahedron illustrated in Fig. [2^ 
gets pushed toward the vertex at — i. The tetrahedron 
in Fig. [2)d, on the other hand, has its upper and lower 
edges pushed together, toward the xy-plane. Both these 
deformations move the spin roots (see [36[) away from 
the north and south poles, which increases the probabil- 
ity that F^ = ±2 in the corresponding spinors, as favored 
by the quadratic Zeeman effect. Now these two types of 
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deformed tetrahedra have different spin-dependent ener- 
gies. In the first case, the magnetization m of the spinor 
becomes nonzero, increasing the /3 term of Eq. ([3|). In the 
second case, the magnetization of the spinor is stiU zero, 
by symmetry, but one can check that the spinor is no 
longer orthogonal to its time-reversal, so 9 is nonzero [55l|. 
Therefore the orientation of the ground state tetrahedron 
is determined by whether the m? term or the \0\'^ has a 
larger coefficient in the Hamiltonian. If the coefficient of 
the magnetization term is very large, then the state with- 
out any magnetization is preferred. The detailed calcu- 
lation shows that the relevant comparison is between c/? 
and 4/3. 

The main source of anisotropy is different at suffi- 
ciently low magnetic fields 26]; the cubic Zeeman effect, 
proportional to B^ , then dominates over the effective 
potential in Eq. ([4]), proportional to B^. Nevertheless 
the the order B'* effect we have calculated can dominate 
over the B^ effect, even when the magnetic field is small 
enough for the perturbation theory just described to ap- 
ply. This is possible because the denominator no/3 in Eq. 
(|¥|) is small compared to the hyperfine energy splitting 
Ahf ■ For spin 2 atoms the effect of the magnetic field is 
given by 



Vm = ^t ^[i^^BY^A^^^^A^^^i^^^^ij (40) 
where Ahf is the hyperfine coupling. It follows that the 

2 d2 

quadratic and cubic Zeeman effects are qF'^ — ^g^ F^ 

and ^g .2 Fz- The analysis we have given applies when 
the magnetic field is weak enough that the wave function 
is not drastically distorted, Eq. ([55]) . but strong enough 
for the second order effect of the quadratic Zeeman term 
to dominate over the cubic Zeeman term. These condi- 
tions combine into 



no 



(3 « fisB « y/AHFnof3 



(41) 



For ^■^Rb at density 5 x lO^^/m^, no/3 = 3 nK and Ahf = 
160 mK while fiB — 67/iK/Gauss; hence the anisotropy 
potential used here is actually valid for a wide range of 
magnetic fields, between .04 mG and .3 G. 



B. Kinetic Energy 

Now that we have estimated the energy due to mis- 
alignments with the quadratic Zeeman field, let us de- 
termine the kinetic energy of vortices in terms of the 
rotations and rephasings, Eq. (|16p . This energy will be 
the source of the repulsion that keeps the vortices apart 
within the vortex molecules. The kinetic energy is de- 
fined by the gradient terms in the Hamiltonian Ti; the 
analogue for liquid crystals is the elastic energy favoring 
alignment of the order parameter. 

In order for the kinetic energy of a vortex to be min- 
imal, the field on a circle around it should trace out 
a geodesic, as mentioned above. For the cyclic phase. 



geodesies take the form given by Eq. ([18)) . To see that 
vortices are geodesies as functions of (j), suppose the field 
of a vortex is given far away by the radius independent 
expression 



^(r,0) = V^F(0), 



(42) 



for an appropriate spinor function F{(j}). The main con- 
tribution to the energy of this vortex is from the kinetic 
energy, which can be estimated by integrating from the 

core radius ac '^ h 

healing length (see 
energy is 



^o . This is the same as the spin 
16j for the definition). The kinetic 



E 



(fr Vi/)'!' • Vt/i 

2m 

""i^'(^)i^^^r2^ 

^'"°ln^r|F'(^)|^d0 
ac Jo 



••dr 



2m 



(43) 



Now the curve parameterized by F{(j)) adjusts itself so 
as to minimize the last integral, while maintaining the 
topology of the circuit traced out by F{(p) in the order 
parameter space. One can show that an integral of this 
form is minimized when F((f>) traces out a closed geodesic 
in the ground state space. The length of a closed curve is 
defined by J^^ \F'{(j))\d(j), so it is not surprising that the 
geodesic of charge F, which minimizes this expression, 
also minimizes Eq. (j43|) . Furthermore, if the geodesic 

has length Ir, then J^'' d(l)\F'{(P)\^ = §^. 

A geodesic in the cyclic order parameter space (which 
has the local geometry of a perfect sphere in four di- 
mensions) can be described by a rotation at a fixed rate 
around a single axis as in Eq. (fTS]). (For a shape less 
isotropic than a tetrahedron, the order parameter would 
rotate around a wobbling axis according to the rigid- 
rotation equations.) Substituting the symmetrical i^((/') 
into Eq. (gS]) gives J4i| 



E : 



h hq , R 

'■ in — 

2m Oc 



dcfii/r {—a- 



27r 27r' 



2^ 



1 ^ r 2/ ^ SI'S + l)x 

■ In — a (n,njQ,j + ) 

47rm flc 3 



(44) 



where the general expression for any spin and phase has 
been given in terms of the quantum fluctuation matrix 
Qij — ;f-(V''^ ' \ ^ ' V') ~ ''3" % and the magnetiza- 
tion per particle Mi = ^ii. The spin 2 tetrahedron state 
appears to be isotropic as long as one does not go be- 
yond second order correlators, as seen from the following 
calculations: 



<Xz\F^F,\xz>^25,, 



(45) 
(46) 



12 



and hence M = Q = Q. Eq. |44] implies that the en- 
ergy is proportional to ^p ~ {6^ + 2a^), a generalization 
of the Pythagorean theorem showing how to combine 
the amount of rephasing and rotation to get the total 
geodesic length. A spin rotation costs twice as much en- 
ergy as a rephasing by the same angle. 

In order to study vortex stability and Coulomb forces, 
let us define the "energy index," 



r \2 



(s> 



(s>^+^(i^>^^ 



(47) 



which is a fraction for each tetrahedral charge from Sec. 
II Bl The energy of a vortex is given as 22LJM j^ il/^(r)^ a 
multiple of the energy of an ordinary phase vortex. The 
force between a pair of vortices can be expressed very 
simply in terms of Ie ■ 

The force follows from an estimate of the e nerg y of 
a cluster of vortices. Using ideas from Ref. [40|, we 
think of the cluster as forming the "core" of a bigger 
vortex, as illustrated in Fig. [HI (We are not necessarily 
assuming that the vortices are bound together.) Draw 
a circle of radius X just around the group of vortices. 
The kinetic energy can then be found as the sum of the 
energies outside and inside of X\ for the case illustrated 
in the figure this energy is approximately 

7r^([Js(X)lnA] + [(/^(l) + /^(2) + /^(3))ln:^]), 

TO Lx CLc 

m 

where Lx is the diameter of the group being combined 
together and R is the radius of whole system. The first 
term describes the energy outside of X. Sufficiently far 
outside of X, the field should have the form of a rota- 
tionally symmetric vortex. The rotation and rephasing 
of this vortex, measured along a circle outside X, is ob- 
tained by multiplying the group elements which describe 
all the individual vortices according to the rule for vortex 
unification Eq. ((22)) . Hence the energy outside X is given 
by an expression like in Eq. (|1^ , except that the integral 
must start at the radius Lx of the circle, so Oc must be 
replaced by Lx- The energy inside X is approximated 
by adding the energies of the three vortices in it, which 
are calculated like in Eq. (|^H) but now with R replaced 
by Lx- 

This approximation makes a small error (compared to 
In ^) by ignoring the region where the vortices' fields 
overlap. This error has a special scaling form if g = 
and if there are no other vortices in the condensate. In 
this case, the kinetic energy gives a complete description 
of the vortices outside their cores, and has a symmetry 
under rescaling. Therefore (as in |46|). one can show that 
the difference between Eq. (P5|) and the actual energy has 
the form 

AE=^-^fi^,^) + Oi^) (49) 

TO Li3 Li3 Lx 

where Ly refer to the sidelengths of the triangle formed 
by the vortices. This correction can be ignored relative 



X 

.2 



X 



3' 




FIG. 6: A group of vortices is combined into the core of a com- 
posite vortex X. Tlie energy outside the blacked-out "core" 
is calculated from the resultant winding number of the vor- 
tices inside, and the energy inside the "core" is calculated 
by looking inside the core to see which vortices are actually 
there. 



to the logarithmically divergent terms kept in Eq. 

as long as the side-lengths all have the same order of 

magnitude, since in this case / has no singularity. 

The general expression for the energy of a set of vor- 
tices with charges F^, after rearranging the formula to 
emphasize the dependence on the diameter of the set L, 
is: 

E^ « ^!^[V/s(FO - /^(]TF,)]ln- 

i i 

+ ^^/.(nr.)ln^ (50) 

i 

As long as the vortex cores are well-separated, the error 
in this estimate depends only on the ratios of the dis- 
tances between the vortices, as in the three-vortex case, 
Eq. (I49p . The energy function has a similar form to 
the two-dimensional Coulomb interaction between vor- 
tices in an ordinary single-component condensate, and 
shows that the vortices repel or attract each other ac- 
cording as the index of the combined vortex is greater or 
less than the sum of the separate vortices' indices. An 
interesting difference is that the interaction is not a sum 
of interactions of all pairs (consider three {A, 0) vortices, 
for example, or almost any other set of three vortices). 

When there are just two vortices, the force between 
them can be calculated by differentiating Eq. ([50]) . giv- 

ing -%a[/^(Fi) + /b(F2) - /^(FiFs)]. This is the 
exact (if L ^ Oc and q = 0) expression for the force be- 
tween the two vortices because the error (see Eq. (09])) 
reduces for two vortices to a constant. Improving on Eq. 
()50|) depends on finding the spin texture in the overlap 
regions, and these equations are nonlinear when not all 
the vortices use the same symmetry axes. 

In the next section, we will apply Eq. (|50p more 
broadly. When there is a small cluster of vortices of size 
L, but there are some other vortices besides, the energy 
of just the vortices in the cluster can be estimated by 
replacing R in Eq. (|5D|) by the distance D to the nearest 
vortex not in the cluster. Likewise, for q 7^ the energy 
of a cluster with a net charge that is field-aligned can be 
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estimated by replacing R by Lq. As long as L <C ig, 
the detailed expression for the energy has the form given 
by Eq. (|49p. but when L approaches Lq, the anisotropy 
energy starts to compete with the kinetic energy and to 
distort the spin texture around the vortices. 



III. 



CHEMISTRY OF VORTICES 



In this section we will first discuss stability of isolated 
tetrahedral vortices and then determine when these vor- 
tices can combine to form molecules (and what the spin 
texture around a molecule looks like). Some of these 
molecules are only metastable and can each break up in 
several ways. 

Bound states of vortices will be formed out of stable 
tetrahedral vortices. These are the vortices based on 120° 
symmetries of a tetrahedron, accompanied by a phase 



shift of ^ or 



-^, and vortices based on 180° symmetries 



without a phase shift. Vortices with larger rotations or 
phase shifts will not occur as components of molecules. 
The net charge of the bound state must be field-aligned, 
with a schematic relation 



Q^l[r, 



(51) 



Thus, if J2ii^'i + 2a?) < (^^ + 2a^), the energy decreases 
as the vortices move apart, so the fragments will move 
apart by themselves the rest of the way. On the other 
hand, the vortex is absolutely stable if 



B|)"-^(£)'>(^)"-^s 



for every set of Ti = (e 



J8 



) such that 1 f Ti = r 

(53) 



Any tetrahedral vortex not satisfying this "absolute 
stability criterion" can break up into a lower-energy 
state, and we will assume that this break-up happens 
spontaneously for these component vortices. For ex- 
ample, the vortex {id, Att) should break up into two 
(id, 27r)'s, halving the energy. (In fact, these two singly 
quantized vortex can break up even further.) Generaliz- 
ing this example, any vortex with charge {g, 9) whose cir- 
culation is bigger than 27r can break up into {g, — 27r) 
and (0,27r) since 9^ > (2^)^ + {0 ~ 27r)2. This leaves 
only finitely many vortices that have the possibility of 
being stable: all the ones with phase winding numbers 
not more than 27r. Some of the vortices with \0\ < 27r are 
also unstable; trial and error finds decay processes such 
as: 



between the field-aligned charge and the component 
charges. A composite vortex is stable if Coulomb repul- 
sions can prevent the collapse of any subset of its com- 
ponents; since the Coulomb interaction is not a sum of 
pairwise interactions, it is not enough that every pair of 
vortices repel each other. Finally any bound state whose 
charge Q — {a, 0)^ has big enough rotation or rephasing 
angles can in principle decay into molecules with less en- 
ergy, but this can only occur if thermal energy overcomes 
the Coulomb repulsion between the component vortices. 
The following expands on this general picture and points 
out some interesting sidelights. 



A. Stable Tetrahedral Vortices 

Only stable vortices will be found in the core-region 
of composite vortices. Since a weak anisotropy term can 
be neglected near component vortices, we can enumerate 
the stable vortex types assuming that q — Q. 

Absolute stability of a vortex with charge F implies 
that if F is forced to break up into the fragments F^, 
then the energy of the fragments grows as they separate 
from one another. The energy of the fragments can be 
found by applying Eq. ^. Note that niTj = F by 
conservation of charge. If F has phase and rotation angles 
and a then 



E 



fragments ~ ^^^V/^^ 



9?+ 2a?) 



9^-f2a2)lln — Vcnst. 



1. {-id, 0) -^ (C, 0) * (C, 0) or {R, 27r/3) * {R, 27r/3) * 
(i?,-47r/3). 

2. {-id, 27r) ^ {R, f ) * {R, f ) * {R, f ) 

3. (i?2, 47r/3) can break up into {R, 27r/3) * {R, 27r/3) 
and also (S'-\47r/3) * (C,0) 

4. (i?2,-27r/3) might break up into (5"\-^) * 
(CO). 

5. (C, 27r) could break up into {R, ^-) * (S'"\ ^). 

All of these decays lower the energy index. There are 
also two vortices which can break up without the energy 
indices changing: 

1. (0, 27r) has the same charge as {R,^-) * {R-^ ,^). 

2. {R,—^) has the same charge as (P^^,— ^) * 



(Q- 



-)■ 



(52) 



The energy may either increase or decrease after one of 
these processes, the logarithmic term which dominates 
the energy (see Eq. ([50])) does not change, so the remain- 
der term needs to be calculated to determine whether 
these break-ups raise or lower the energy. A point vortex 
whose charge is (0, 2n) is probably unstable if a ^ /3 
because a pure phase vortex has to have an empty core 
(with energy density of order an^), while the two frag- 
ments it breaks up into just have non-cyclic cores (energy 
of order (in^). [Sq 

We can now enumerate the vortices which are stable, 
noting that some vortices (such as (P^, ^)) can break up 
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similarly to the ones just listed on account of symmetry. 
The only vortex types that have not been eliminated are 
the vortex with one-third circulation, {R, 27r/3) (energy 
index i) and the currentless vortex (C, 0) (energy index 
i) as well as possibly (i?, — ^) (energy index |), and 
also the inverses and conjugates of these. 



B. Vortex Molecules at g 7^ and their Spin 
Textures 

In this section we will describe a qualitative wave func- 
tion for the example in Section II Al to illustrate how the 
field of a vortex molecule deforms in response to the 
anisotropy energy as one leaves the region containing the 
two vortices. We will then give the binding criteria which 
describe how to use group theory to check which sets of 
the stable tetrahedral vortices (from the previous sec- 
tion) can combine together to form a vortex molecule at 
a nonzero magnetic field. 

A vortex with a charge that is not compatible with the 
magnetic field has an energy that grows proportionally 
to the area of the condensate. Eq. (P5)) shows that the 
vortices in a cluster can avoid this energy cost if they have 
a net charge corresponding to a rotation around the z 
axis. More specifically, a set of tetrahedral vortices, with 
topologies Ti , has to have a net topology of the form 



nr. = (i?™,^ + 2.n). 



(54) 



The "molecule" can then have the rotational and phase 
windings [a, 6)3 = (^ + 47rj,27r(f + n))3 at infin- 
ity, for any j. One example of a molecule for the Mqs 
phase is the composite vortex discussed at the beginning, 
{A, 0){A, 0), which partner up to make a compound vor- 
tex with Q = {2tt, 0)3 in the field-aligned condensate. 

Writing a wave function that describes this example 
even qualitatively is a little more complicated than it 
sounded in Sec. II Al A possible wave function is illus- 
trated in Fig. [71 We can build this wave function up 
in stages. The most obvious attempt at writing a wave 
function fails to eliminate the quadratic energy cost: 



V': 



'misalign 



,-f(0i+02)(VfF, + ;ijF-) 



"0X3, 



(55) 



where (f)i — arctan ^_^' is the polar angle measured with 



,-th 



respect to the location {xi,yi) = ±(2-,0) of the i"" vor- 
tex. The two 180° rotations about A combine into a 2tt 
rotation on a circle surrounding both vortices, but not 
about the magnetic field axis. This wave function takes 



^^^iV2F^ 



'F,: 



noxs at infinity where 



the form ip ^ e 
(j)i and (j)2 approach (j). The tetrahedra rotate around 
a tilted axis, so they do not stay aligned appropriately 
with the magnetic field except on the x-axis and the y- 
axis (where the tetrahedra are reversed, but this is still 
a ground-state). The Zeeman energy of tpmisaiign is still 
infinite. This wave function agrees pretty well with the 



one illustrated in the figure within the region surrounded 
by the large circle. The vortex cores are cordoned off 
by the small circles, where the 180° symmetry axes are 
marked by dots. Note that the tetrahedra three rows in- 
side the large circle behave like the tetrahedra at infinity 
described by Eq. (|55)) : they rotate through 360° about 
the order-two axis bisecting their right edge, producing 
all sorts of arbitrary orientations. 

Another attempt, which uses the R axis instead of the 
A axis to eliminate the tilting of the tetrahedra. 



Ydiscontinuous — ^ 



"0X3, 



(56) 



is a complete fiasco, since this function is not continuous 
along the line connecting the two cores. (The R axis has 
the wrong symmetry, so as (x, y) circles around {x2,y2)^ 
i^ discontinuous chaugcs from X3 to e'^-'^^'Xi 7^ X3-) 

Luckily, since a 2-k rotation around one axis can be 
deformed to any other {2-k rotations all correspond to 
—id G SU2) one can produce a spin texture which has a 
finite Zeeman energy and is continuous. A hybrid of Eqs. 
([55)1 and (|56p which achieves this, illustrated in Fig. [3 



V' 



-#(_F, sin M(r)+_F, cos Kr)) pi't'^^^^^^f^ 



■i{'t>l+<t>2)U 



\F. + ^F^) 



«0X3 



where ^{r) must satisfy 

/x(oo) = 
For example, we could define 



arccos \ / - 



cos/i(r) 



D^ 



3D2 



(57) 

(58) 
(59) 

(60) 



where _D is a variational parameter. (Probably D ^ L is 
optimal.) 

This field has two vortices at (± -j , 0) . The texture 
varies continuously, except at the at the very centers of 
the vortices. (In the figure, the cores are inside the small 
circles. The tetrahedra's orientations change rapidly 
around the circles' centers, but each tetrahedron inside 
has a similar orientation to the tetrahedron next to it out- 
side the circle.) Around the cores, the tetrahedra rotate 
around a tilted h' axis which is aligned with the symme- 
try axis of the tilted tetrahedra around the vortex; the 
tilting is carried out by the rotations described by the 
first two factors in Eq. ([57)1 . just as discussed in Section 
II Al Near vortex 1, the wave function is approximately 



■0 « e 
and near vortex 2, 



~{^F^+^F^: 



«0X3 



(61) 



0: 



-i7r(F^sin/i(i)+F^cos^(i)) 



F, + ^/i"F,) 



"0X3- (62) 
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(The local symmetry axis ft' is the rotation of the 



standard axis 



3^ 



y, through a half-turn about 



^sin/i(-|) + ycos/i(-j), according to Eq. (fTB]) .) 

Though the two tetrahedra are rotated into misaligned 
positions inside the large circle, the excursions from M.q3 
to M. are brief and therefore cost only a finite amount 
of energy. Indeed, the first two factors fix the field up at 
infinity by applying a continuously varying rotation to 
the overall texture. These factors change the axis from 
A to i? at large r as one can see by replacing 0i « (f>2 by 
(f) and using Eq. ([55)1 : 



^ 



■.-i4>F, 



"0X3- 



(63) 



The 360° axis changes from the A axis to the R axis 
as one crosses through the transition region indicated by 
the large circle in Fig. [T] As you follow a radius out- 
ward past the circle, the tetrahedra are tipped by dif- 
ferent amounts. (The tetrahedra on the negative x-axis 
have to be tipped the most, though their original ori- 
entation is compatible with the magnetic field! The face 
which is on top is changed. This reorientation is required 
to make the amount of tipping continuous: the tipping 
angle increases more and more as goes from to tt.) 

The first two factors manage to fix the orientation at 
infinity without introducing discontinuities like in Eq. 
([56)1 . Although they seem to have a vortex-like discon- 
tinuity at the origin, their discontinuous parts cancel on 
account of Eq. (|58|1 . (Hence the tetrahedra near the ori- 
gin in Fig. [7| all have roughly the same orientation.) In 
short, though (^,0)^ = (— id, 0), the cancellation of the 
Zeeman energy at infinity is not automatic! The topolog- 
ical classification just implies that the field in Eq. ([55]) 
can be deformed as r ^ oo so that the Zeeman energy is 
small. 

The kinetic energy of this composite vortex, according 
to Eq. dnni), is 



(64) 



m ac 

■nh^TiQ ^ ,, . , ^.., R 
+ -lE{{-id,Q))\n—. 

m ttr 



Since Ie{{A,0)) = i and /£;(— id, 0) = 2, the coefficient 
of the first term is negative, so the vortices are driven 
apart in order to decrease the kinetic energy. This effect 
competes with the anisotropy energy, which tries to bring 
the vortices together. The net energy is therefore 



E « fiD/L) 



2n2 



q'D 



.^!^(21n:^+ln^). (65) 
m L a 



This is a more complete version of Eq. ([5]) . / summarizes 
the dependence of the anisotropy energy Eq. ^ on D 
and L, but is too complicated for us to figure out! 

The competition between the attractive Zeeman-term 
and the repulsive kinetic-energy term determines the size 



of the vortex. Assuming D = L and minimizing over L 
gives 



ipK^P 



(66) 



This composite core-size is also the scale for the de- 
cay of tipping of tetrahedra due to the competition of 
anisotropy and kinetic energy (see Sec. Ill A|l . just as the 
core size of a spin vortex in a spin 1 condensate is equal 
to the magnetic healing length. Vortex molecules will 
have a size of the same order of magnitude (even if there 
are more than 2 subvortices) except when the component 
vortices have a short-range repulsion, like the two exam- 
ples at the end of Section fill Al (Such components form 
smaller molecules; vortices never form larger molecules. 
If a pair of tetrahedral vortices is stretched beyond Lg, a 
"cord" of tipped-tetrahedra forms between them, a sim- 
ple version of the string imagined to hold the quarks to- 
gether in a rapidly-rotating baryon.) 

To find the total energy of the molecule, substitute L 

into Eq. (|65p and simplify using ac = 



4-R . 
nQf3m ' 



_E,, 



27r— — Ini?- 



(67) 



We have dropped the contribution from the quadratic 
Zeeman term since it adds something that is independent 
of q. 

This formula can be compared to the energy one ex- 
pects for a simple vortex with the same charge (27r,0)3: 



Evortex = lEi{'2n,0)3)n In — 

m a' 



(68) 



which includes a core energy ec and allows for a physical 
definition of the core size. If we take a' to be the size of 



the molecule, Lg ~ - 
to be 



"°^ ' , then the core energy has 

(69) 



TT in 

m q 



The core energy becomes large as g -^ because the 
kinetic energy of the component vortices diverges as the 
logarithm of the core size Lq. 

It makes sense to regard this composite vortex as a 
molecule because Coulomb repulsion keeps the vortices 
in the core separate. The previous discussion can be gen- 
eralized by listing a set of binding criteria; these ensure 
that a set of tetrahedral vortices will form a stable or 
metastable composite vortex: 

1. Each component vortex is one of the stable q ~ 
vortices from Section fill Al 

2. The kinetic energy is not decreased when any sub- 
set of the component vortices coalesces into a single 
vortex. 
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FIG. 7: An illustration of the composite vortex {A, 0)*{A, 0). The magnetic field is perpendicular to the figure, so the tetrahedra 
prefer to be oriented with a face or a vertex facing up (since c > 4). The two small circles enclose the 180° vortices, with 
symmetry axes indicated by black dots. The large circle indicates the transition region where the axis of the 360° rotation 
changes relative to the tetrahedra, from the A to the 7? axis. The figure uses a functional form with a rapid jump for p(r) 
(unlike in Eq. (|60|) ) for simplicity, so that the transition region does not overlap the cores. For the probably more realistic 
form given by Eq. (|60p , the tetrahedra are already tipped at the centers of the vortices so that the local rotation axes n^ do 
not have the standard orientations illustrated in Fig. [2^. 



3. There is no way for the component vortices 
to form submolecules that can break apart. 
This would occur if the components could be 
rearranged and then partitioned into r sets 



{ri,r2, . . . ,rji}, {rji+i,rji+2, . . . ,rji + j2}, ■ ■ ■ , 

such that each subset forms a molecule that 
is compatible with the magnetic field (i.e., 
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(i?^,27r(^ + Uk))) and 



njlH hifc p 

such that t/ie sum of the energy indices of these 
suhmolecules is less than the energy index of the 
original molecule. 

The vortex molecule {A, 0) * (A, 0) clearly satisfies all 
these conditions: Condition 1 is satisfied because (A,0) 
is one of the stable vortices found in Section IIII Al Con- 
dition 2 is satisfied because the vortices repel each other. 
Condition 3 is easy to check for a diatomic molecule 
like this one, since it can only break up into individual 
"atoms" ; neither of the fragments [A, 0) is compatible 
with the magnetic field. 



C. Metastable Vortices and How They Decay 

Not all of the vortex molecules satisfying the three con- 
ditions above are absolutely stable. The analogue of the 
absolute stability condition, Eq. ([55]) . also selects a fi- 
nite set of aligned vortex types when g ^ 0, this time 
from among the group elements listed in Eq. (|30p . The 
absolutely stable charges are 



27r 2^ 



±( 



,27r 47r 



(70) 



For any other vortex topology Q — (q, 0)3, one can find 
vortex topologies Qi such that Y\^ Qi — Q and 



IE{Q)>J2^E{Q^ 



(71) 



Point vortices with such a topology Q would likely break 
apart spontaneously. 

There is actually another pair of charges that could be 
absolutely stable, but the energy index estimate is not 
accurate enough to decide the issue: 



(0,±27r)3,(±f,Tf)3 



(72) 



These vortices can break up into pairs without changing 
the net energy indices, reprising the ambiguous behavior 
of the two vortices at the end of Sec. IIII Al (See Section 
lIVBI for an answer.) 

There are some composite vortices of other charges be- 
sides the four listed in Eq. (|70p which are long-lived. The 
absolute stability criterion misses this possibility because 
it ignores the details of the vortex cores, drawing all its 
conclusions from the topology of the vortices far away: 
Suppose the initial vortex Q is a cluster of vortices with 
topologies ri,r2, . . . ,rr (see Fig. 18]). The decay prod- 
ucts discussed in the previous paragraph, Qi, may be the 
combined topology of other vortex clusters. The energy 
of the vortices after the reaction (Fig. [5Jd) is smaller if 



Einit = Ie{Q)it In — 

m L 



> 



\^ 



[IE{Q^)^T 



n^h^ 
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= E 
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(73) 
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FIG. 8: Illustration of the absolute stability criterion. A 
composite vortex (a) and a possible set of composite vor- 
tices (b) it can break up into. Even if Q = Q1Q2Q3 and 
Ie{Q) > Ie{Qi) + Ie{Q2) + lEiQs) so that the energy would 
decrease, the break-up might not occur spontaneously. If the 
component vortices in (b) are different from the component 
vortices in (a), then the vortices making up Qi, Q2, and Q3 
would have to be produced in a "chemical" reaction from the 
components of Q. 



The terms proportional to In j- stand for the kinetic en- 
ergies outside the composite cores, and the terms propor- 
tional to In — stand for the energies within the compos- 
ite cores; the latter contributions do not matter once the 
composite vortices are far apart (i? 3> L). Therefore, the 
energy-index relation, Eq. (|7ip . implies that the energy 
decreases when Q ^ Qi* Q2* ■ ■ -Qr- 

However, at zero temperature, a vortex molecule sat- 
isfying the binding criteria cannot break up, even if the 
total energy would end up smaller. Although the meta- 
morphosis of Fig. [8^ into Fig. ^p lowers the energy, the 
process will not occur spontaneously. According to con- 
dition 3, the vortices in Fig. [Hja are different than those 
in Fig. [8^. And according to conditions 1 and 2, there 
arc no spontaneous chemical reactions that can occur to 
make the components in Fig. [5Jd out of those in Fig. [5^. 
Thus, at zero temperature, composite vortices besides 
the ones with charges listed in Eq. (jTO]) . whose energy 
indices seem to be too big, can still be stable. At nonzero 
temperature, such a molecule will only be metastable be- 
cause it can decay after reactions inside its core produce 
the vortex types that appear in Fig. [8)3. These reactions 
are prevented by energy barriers, so the decay will occur 
only after a long time. 

The vortex molecule made up of {A, 0) * {A, 0) is an 
example of a metastable vortex; its charge (27r,0)3 did 
not appear in Eq. ([70]) because it is the same as the net 
charge of the three field-aligned point vortices (^, ^)3 * 
(¥' ~T^)3 * (¥' ¥)3- ^^^ energy index of the molecule 
is 2(1^) -I- 0^ (see Eq. (liTt ) while the energy index of the 
three point vortices is smaller, ^ + | + |- Nevertheless, 
since these three vortices are not present in the core of 
the original vortex the decay cannot occur spontaneously. 

A vortex with a composite core can break up by a 
combination of vortex-fusion and vortex fragmentation. 
Some subsets of the original component vortices fuse and 
the fused vortices each break up into some other vortices. 
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These regroup into clusters each of which has a charge 
compatible with the magnetic field. Then each cluster 
goes its own way. (The fusion and fragmentation steps 
might sometimes happen more than once.) Conditions 
1 and 2 ensure that at least one of these steps will be 
opposed by the Coulomb potential, but the total energy 
will decrease if the energy index decreases. The vortex 
molecule will be long-lived because its components do 
not know that the hard effort of fusing will allow them 
to change into vortices which can separate. 

Trial and error yields a couple of ways in which the 
{A, 0) * {A, 0) bound state can break up. One possibility 
begins with the two component vortices coalescing, 

Oqr Oqr Ajr 

{A, 0) * (A, 0) ^ i-id, 0) ^ {R, —) * (i?, -) * {R, - — ) 

(74) 
The other begins when one of the components breaks up. 



27r, 



27r, 



{A, 0) * {A, 0) ^ [iR-\ - — ) * (g, — )] * {A, 0) 

-(i?-\-|^)*[(g,^)*(A,0)]. (75) 

In the first process, the two vortices come together, in- 
creasing the kinetic energy in accordance with Condition 
2 (as shown by calculating the /b's and substituting into 
Eq. (j50[) ). The resulting vortex breaks up into three vor- 
tices which can separate from each other because they 
are compatible with the magnetic field. The increase in 
energy during the first stage is given by Egi — Einu — 
Trlio^ In h. where E^i = -K^^^lEi-id, 0) In -^ is the en- 
ergy of the intermediate vortex. Thermal fluctuations 
have a chance of driving the (yl, 0)'s together, in spite of 
the energy increase Ea\ — Emit- 

In the second process, Eq. ([75]) . the (A, 0) vortex first 
splits up into two vortices. The first of these, {R~^, —^) 
is compatible with the Zeeman field and can leave. The 
remaining two vortices form a new molecule which can- 
not break up because Q is a rotation around the wrong 
order 3 axis. For this process Condition 1 requires that 
the energy increases during the initial fragmentation. To 
check this, note that the energy of the intermediate state 
is 



E, 



s2 






lEHd,0)lnj} 



lEiQ,^)+lEiA,0)]lli- 

6 ac 

(76) 



and the energy barrier is 



E 



s2 



:n^[{lE{R-\-^)+lE{Q,^)-lE{AM^^- 



m 



irnQh^m L 

T In-. 

D a 



(77) 



This energy barrier is lower than Esi — Ei^it , so Eq. (I75p 
is a more common break-up route. (In a finite conden- 
sate, thermally excited break-ups can be observed only 



if a vortex molecule is somehow prevented from wander- 
ing to the boundary of the condensate and annihilating 
before it can decay.) 

These two examples illustrate the meaning of the sta- 
bility conditions. Conditions 1 and 2 ensure that frag- 
mentation and fusion processes cannot happen sponta- 
neously. The third condition simply points out that 
vortex clusters like (A, 0)^ * {B,0)^ will not be stable 
because the components can sort themselves into field- 
aligned groups and break up without any thermal assis- 
tance. The second condition can be difficult to check 
for a composite vortex with three or more sub-vortices. 
One must consider subsets of every size and check that 
they cannot lower their energy by collapsing all at once 
into one vortex. Just knowing that any two vortices of 
the subset repel each other does not guarantee that the 
set of vortices do not collectively attract each other! An 
example is the set of three vortices {A,0)^. Any two of 
these vortices would combine to form a vortex with a 
360° rotation, {—id, 0), whose index 2 is higher than the 
sum, 5 + 5, of the indices of the collapsing pair. On the 
other hand all three vortices could form a vortex {A~^,0) 
with energy index i < 3/_e(^, 0), so the three vortices 
can collapse simultaneously to lower the energy. An even 
more counterintuitive complication is that, because of 
the noncommutative behavior of the combination rules, 
more complicated fusion processes can occur. A vortex 
can change its type by circling around one vortex so that 
it can fuse with another vortex. (See Appendix IB] ) For 
a bound state of many vortices there will be many possi- 
bilities for how the vortices meander around each other 
before some of them fuse. To test Condition 3, one also 
has to enumerate all possible wanderings. 

A mathematical statement of our results is that there 
will be (at least) one solution to the time-independent 
Gross-Pitaevskii equation for any set of vortex topolo- 
gies that satisfy the three conditions. If q is small, the 
energy at the top of the barrier is greater, by a loga- 
rithmically large amount, than the energy of an initial 
variational state like the approximate wave function Eq. 
([57)1 . A solution to the Gross-Pitaevskii equation should 
result if one starts from this qualitative texture and lets 
it relax to a local minimum of the energy. There is not 
enough energy for the wave function to get over the en- 
ergy barrier, so the wave function should get stuck in a 
local minimum. (The energy of the intermediate state is 
not known precisely because of the rough estimates we 
have made of the Zeeman energy and the kinetic energy, 
but these errors are small compared to the height of the 
barrier.) 



IV. ADDITIONAL EXAMPLES 

Now we can construct some other, more interesting, 
examples. We will use the algebra of the group of vortex 
charges to find molecules whose net charge is interesting 
in different ways, and we will use the energy index to test 
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Components 


Net Charge 


c 


Stable? 


1 


(^,0)*(^,0) 


(27r,0)3 


c>4 


Metastable 


2 


(p-\^)*(0-\4f)*(i?-\if) 


(0,47r)2 


c<4 


Metastable 


3 


(Q,4f)*(A0) 


( 3"'— )3 


c>4 


Stable 


4 


None 


(47r, 0)2,3 


Any value 


Metastable 



TABLE I: Examples of vortex molecules. The tetrahedral charges of the components of the molecules and the net aligned 
charge are given. The condition on c determines how the tetrahedra are oriented far from the vortex, due to the magnetic field. 
The final column indicates whether the vortex molecule is expected to have the absolute minimum energy of all vortices with 
a given net aligned charge. The second molecule might actually not be bound-see the text. 



whether they are stable. 

The parameter c will be less than 4 for some of these 
examples. If c < 4 the ground state orientation of the 
tetrahedron will be as in Fig. [^b so the z-axis is an order 
2 axis. The aUgned topologies have the forms (7m, 27rm)2, 
as described in Section [TBI Of these, the only absolutely 
stable topologies are ±(7r,0)2, ±(0, 27r)2 (and possibly 
i±7r,±2nh). 



is satisfied. Condition 3 is clear. To check condition 2, 
let us first consider whether one of the pairs of vortices 
in the trio can coalesce. Using conservation of topolog- 
ical charge helps to avoid enumerating all the ways the 
vortices can braid around each other. If the first two 
vortices have coalesced into a vortex {X, ^) (after some 
permutation) and the third vortex, by winding around 
the other two vortices as they collapsed, has changed to 
(y-\^), then 



A. A Doubly Quantized Pure Phase Vortex 

First let us find a vortex molecule whose phase winds 
by 477. In single component condensates, such vortices 
are usually unstable; one has been observed to break up, 
maybe into an entwined pair of 27r vortices[30J. If phase 
and spin textures were completely independent of one an- 
other, doubly quantized vortices would not be any more 
stable in the cyclic condensates; but fractional circula- 
tions are "bound" to certain spin textures (see Sec. IIB[) . 
// we assume the vortex {R, —3^) is stable (at the end of 
Sec. nil Al we could not decide), then a doubly-quantized 
vortex can occur in a cyclic condensate when c > 4. It 
consists of the three parts 

(P-,f)(Q-,f)(i^-,f). (78) 

The phase changes by 47r while the orientation of the 
tetrahedron does not change at infinity as we can check 
using the coordinate system from Fig. [lb- The three 
group elements are 

P^' = ^{^ + ii^x - ^y - ^z)) 
Q'^ ^ ^{l + i{<Tx + cry + a^)) 

R^^ = 2'^"'" ''' '^^'^^ - Cfx- CTy)) 

and their product is the identity. 

Let us discuss the conditions for binding. We have 
not checked Condition 1; it is not easy to check be- 
cause (P~^, ^) has the same charge and energy index 
as (i?, ^) * (Q,^); an accurate solution for the spin 
texture around this pair of vortices is needed. Besides, 
(P~^, ^) might be stable for some ranges of c values, 
but not others. Let us therefore hope that condition 1 



(XY-^,Att) = {id, An) 



(79) 



by conservation of charge. Hence X = Y. Also, braid- 
ing one vortex between other vortices can only conjugate 
its group element. Therefore, Y, like R, is a counter- 
clockwise rotation through 120°. Since X — Y, the ro- 
tation part of the coalesced vortex (X, ^) also is a 120° 
turn and thus the energy index of this coalesced vortex 
is 2 X (1/3)^ + (4/3)^ = 2, which is greater than the sum 
of the energy indices of the two vortices which formed 
it. Therefore the two vortices cannot coalesce sponta- 
neously. (This argument can be generalized to any trio of 
vortices ri,r2,r3 each of which commutes with the net 
charge F. Fusing two of the vortices gives the same result 
(up to conjugacy) no matter how the vortices are mixed 
around first; so braiding cannot make a repulsive interac- 
tion between two vortices into an attractive one.) Finally, 
the three vortices cannot coalesce simultaneously because 
/£(0,47r) > lEiP-\ f) + Ie{Q-\ f ) + Ie{R-\ f ). 



B. A Vortex Molecule which is Stable 

Returning to the original assumption, c > 4, where the 
ground state orientation is illustrated by Fig. [5^, we can 
show that the second charge in Eq. ([72|) does correspond 
to a completely stable vortex molecule. In fact, consider 



27r 

(Q,y)(AO). 



(80) 



This molecule, one of the decay products in Eq. ((75)) . has 
the topology (i?^^, ^) or (using the notation appropri- 
ate for the field-aligned tetrahedra outside the composite 
core), (— ^, ^)3- The energy of this molecule is approx- 
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imately 

^^[Ie{R-^, — ) In — + Ie{Q, — ) In -^ + Ie{A, 0) In 



Ln 



formula that shows how a 47r rotation-vortex can relax 
in the absence of a magnetic field: 



^{cj)-t) 



— i(f){Fx sin TTt-^-Fz cos irt) —iF~ ^ 



"oXo, 



(82) 



noTrft^ R I Lq, 
(In ftln— . 



(81) 



where Lq is the size of the composite core, given by Eq. 

This molecule answers a question from Section IIII CI 
Are there stable vortices with charge (—3^,3^)3? The 
two vortices (— X'^lf^^ * (^T"' T" ^3' ^^'^'^ the same 
net topology as a vortex of charge (— ^, ^)3, and they 
have the same net energy index. Now we can check that 
the composite vortex {Q, ^){A, 0) is a stable realization 
for the charge (—3^, ^)3- Its energy is lower by a finite 
amount than the energy 



irnQiTh ( 1 



{\ + |)ln^ of the pair 



of vortices. This finite binding energy is 






llii^. 



To take another point of view, the minimum-energy 
spin texture with the topology (—3^, 3^)3 imposed far 
away has an asymmetric structure: it has two "singular- 
ities" with topologies (Q, ^) and (A, 0) at a distance of 
order Lq. By contrast, when the topology imposed at a 
boundary corresponds to unstable vortices, the ground 
state has singularities whose spacing is on the order of 
the size of the system R. E.g., in a scalar one might try 



to impose ^{R, 0) 



noe 



2id 



The spacing of the vor- 



tices in the energy minimizing wave function grows with 
R, reflecting the fact that these vortices would repel each 
other to infinity in an infinite condensate. 



C. A "Bound State" of No Vortices 

The final example shows that point vortices are not 
necessary to hold a core together-there is a "composite" 
vortices without any components! In other words we can 
construct a vortex for which the order parameter stays 
in AA . There is still a "composite core" where the tetra- 
hedra leave Mq and are no longer aligned with the field 
axis. The trick is that the amount of rotation in a tex- 
ture around a vortex is defined only modulo 47r, in the 
absence of a magnetic field (because of the SU2 charge 
classification). When B is turned on, the spin part of 
the order parameter space Mq has the same topology 
as a circle, so each additional winding by 27r changes the 
topological charge. (A texture which rotates by 47r about 
the field axis can relax only by using axes perpendicular 
to the magnetic field.) Thus a 47r-rotation-vortex is sta- 
ble in a magnetic field, but since it has zero tetrahedral 
charge, it does not have to have point vortices inside of it. 
(Another way to say this: Eq. (|54p does not uniquely de- 
termine the aligned topology, because R^ = id. Hence a 
(47r, 0)3-vortex can be made from 6 R vortices (i.e., some 
{R, 3^)'s and (i?, — ^)'s) or out of no vortices at all!) 

Such vortices occur for both the c > 4 and c < 4 cases. 
A variational wave function can be constructed using the 



where xo is an arbitrary cyclic spinor. At each moment 
of time t, the expression describes an r-independent tex- 
ture as a function of (j). When i = 0, there is a vortex 
which is a full rotation through 720°. By the time i = 1, 
this vortex has completely dissipated. When q ^ a, iir 
vortex cannot relax in this way because the tetrahedra ro- 
tate away from the orientation preferred by the magnetic 
field before returning to the preferred orientation at the 
end. But Eq. ([82]) has a reincarnation as the description 
of a {Air, 0)2,3 vortex. We replace the time coordinate by 
a function of the radius to give a spin texture that winds 
through Att at infinity but does not have any singularities 
at 0: 



ip = tp{(l); 



1 



1 



ii-r^' 



(83) 



If Xo = X2 or X3, then this wave function, at large r's, 
has the winding number (477,0)2,3. At small r's, the 
wave function is (/)-independent, giving a continuous and 
"coreless" wave-function. (The exact solution not only 
has a more complicated r-dependence, but also a less- 
symmetrical (/)-dependence.) The region r < L is the 
composite core of this vortex in the sense that -0 G -^ 
rather than Aiq. The optimal size L of this region is 
again Lq, as balancing the kinetic and Zeeman energies 
shows. 

This vortex cannot disappear because the classification 
of vortices at nonzero q implies that a = 47r is conserved. 
Furthermore, though it does not satisfy the absolute sta- 
bility criterion, since two (27r, 0)2,3 's have a smaller en- 
ergy, it is obviously metastable-there are no vortices in 
the core to break apart! The vortex can only break up 
if thermal energy causes a vortex-antivortex pair to nu- 
cleate in the core. Suppose a pair involving rotations 
through 27r in opposite directions appears. These vortices 
initially attract each other but if the thermal fluctuations 
pull them to opposite sides of the core the nonlinear cou- 
pling with the background field switches this force from 
attractive to repulsive, and the vortices can separate the 
rest of the way by themselves. 



CREATING AND OBSERVING VORTEX 
MOLECULES 



Let us discuss the conditions under which the Zeeman- 
effect bound states might be observed and the methods 
one can use for observing them. First of all, we must 
assume that q <C /Stiq in order to justify neglecting q 
near the tetrahedral vortices and to justify the pertur- 
bation theory of Sec. Ill Al This is not just a technical 
assumption: above a certain magnetic field the compo- 
nent vortices probably merge. To estimate the maximum 
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magnetic field note that q is related to the hyperfine split- 
ting Ahf via 



\<1\ = 



a4^ 

8A 



HF 



(84) 



for rubidium and sodium atoms 2], and similar relations 
hold for other atoms. Also note that the spin independent 
interaction is 



Airh^ 



(85) 



where a ~ 50 A is an average of the scattering lengths 
corresponding to different net spins and that the spin 
dependent interaction is 



P = 



Anh^Ac 



m 



(86) 



where Aa ^ 1 A [33, [4l| depends on the differences be- 
tween the scattering lengths [3]. The condition for our 
analysis to be applicable, q << no(3, therefore implies 



B «B 
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Ma 



'Ah.^^^^, (87) 

^B V m 



about . 1 G for a condensate of rubidium atoms with den- 
sity no = 5 X lO^'^/cc. 

In order to observe vortex bound states, one might 
start with a condensate prepared with a spin order other 
than the ground state and then watch it evolve as in Ref. 
[16| . Thermal (and less importantly quantum) noise will 
produce perturbations that grow exponentially, produc- 
ing complicated patterns. If the magnetic field is small 
enough, vortex bound states might be found after some 
time. To test whether these vortex bound states behave 
in the way we have been describing, one would have to 
identify the topological charges of the vortices. One could 
then check that vortex sets whose net charge is compat- 
ible with the magnetic field have a size on the order of 
the theoretical value, Lq. One may have to use statisti- 
cal correlations if too many vortices stay around. (One 
could also take a more deliberate approach, choosing vor- 
tex types and imprinting them as in [43| or [2^ . One can 
then observe the subsequent evolution of the vortices to 
see whether they bind together.) 

In fact, identifying the vortices that appear in a spinor 
condensate is difficult; vortex cores in a spinor conden- 
sate are not empty like the vortices in an ordinary con- 
densate; they have nearly the same density as the rest of 
the condensate 48] . One thus has to measure something 
about the spins to observe the vortices. Two possibili- 
ties have already been invented. One can either measure 
the magnetization field as in [I^ or use Stern-Gerlach 
separation to measure the density of the different spin 
species. 

Measuring the magnetization as a function of position 
is less informative for a nonmagnetic phase like the cyclic 
phase than for the ferromagnetic phase studied in Ref. 



|16j . The magnetization outside the core of a vortex, 
where the spinor state is approximately a rotation of the 
unmagnetized cyclic state will be close to zero (see Eq. 
PS)) ). but inside the core, where the order parameter 
leaves the ground-state space M, the magnetization can 
be nonzero. Measuring the magnetic moment in the core 
of a vortex helps to determine the topological charge of 
the vortex. (The magnetization will not provide any di- 
rect evidence of the rotating orientation of the tetrahe- 
dral order parameter, though.) Any vortex one might 
have to identify involves a rotation about an arbitrary 
axis n' as in Eq. p^ . or Eq. P^ which is more conve- 
nient for understanding what a vortex will look like. The 
latter description starts with a vortex whose rotation axis 
is special-say it is parallel to z, and applies some overall 
rotation to it. 

For example a vortex of type {R, ^) can be obtained 
from a vortex whose axis is n = — z. Eq. (J19p implies 
that the vortex is described by 



ip{r,<l)) = D{R)e3 



i(i+-F'.) 



^(/M\/^,0,0,g(r)y|o 



DiR)V^{ /(r)e*^/-,0,0,5(r)\/-,0 



where /(r) and g{r) are appropriate functions approach- 
ing 1 at infinity and i? is a rotation that moves n to n'. 
(The phase, ^, does not matter.) If n' = n = —z, then 
R is the identity, so mx{r) — rny{r) = and 



m^{r,(j)) = ^{r,4>)'^F^^{r,(f) 

= \[f{rf^g{rf]n,. 



(89) 



The magnetization is parallel to the symmetry axis and is 
given by m= |no[5(r)^ — /(r)^]n. Applying an arbitrary 
reorientation R changes the magnetization axis and the 
symmetry axis in the same way, so the general result is 



m= -^Mgirf ~ f{r, 



^W. 



(90) 



Far from the core, the magnetization vanishes. Inside the 
core, the magnetization can be found by noting that the 
top component of the vortex Eq. ((88l) must vanish at 
r = in order to be continuous: 



/(O) - 0. 



(91) 



Since a is much larger than (3 and 7, the density of atoms 
will be almost uniform across the whole vortex and hence 
\f{fY + |5('')^ ~ 1- Eq. (Ell) therefore implies 



m 



(92) 



Hence the magnetization, Eq. (j90p . is approximately 
non' in the core; the atoms have a single unit of hypefine 
spin in the direction of the vector from the center to the 
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FIG. 9: Densities in the five spinor components around an 
order 3 vortex, with a randomly oriented local axis n'. On 
the left is the pattern one might observe experimentally; the 
lighter regions correspond to regions with fewer atoms. On 
the right are plotted the percentage of atoms for each value of 
Fz at some fixed distance from the vortex core. The phases 
and amplitudes of these oscillations should help to determine 
the direction of the local axis. 



fixed vertex of the rotating tetrahedra near vortex the 
core [57]. The inverse vortex, {R~^,—^), has the same 
magnetization (it does not change sign). On the other 
hand, similar arguments show that {R, — %-) will have a 
magnetization approximately equal to —2noh at the core 
center because it is the m — —1 component of the spinor 
in the analogue of Eq. (|88|) that has the phase wind- 
ing for this case. The third stable vortex type, (A, 0), 
will not have any magnetization in its center and would 
be hard to see using this method. Measuring the mag- 
netization reveals vortices of order three but does not 
distinguish between vortices and antivortices and does 
not even indicate the presence of an order two vortex. 
(One can observe the composite vortex described in Sec- 
tion |TVX1) Another deficiency is that the cores are only 
about 1 /im across, so the vortices might be hard to ob- 
serve directly by this method. However, one could first 
allow the condensate to expand in the transverse direc- 
tion so that the atomic interactions decrease. The vortex 
cores would expand; as in experimental observations of 
vortices in single-component condensates, the depleted 
region in / or 5 (whichever corresponds to the compo- 
nent of the transformed spinor with the phase winding) 
would fly apart and the magnetized core would become 
much larger. A magnetized ring would form at the edge 
of the core where the atoms of one magnetization accu- 
mulate more than the atoms of the other. 

The Stern-Gerlach method gives more information 
about the vortices. Though the density is depleted at the 
center, the field around a vortex in a single-component 
condensate is not observable, unless one reconstructs the 
phase variation of the condensate, perhaps using the 
technique described in Ref. [4§|. But in a spinor con- 
densate, the spin vortices produce observable patterns 
in the condensates' Stern-Gerlach images. These images 
capture separately the density of atoms in each of the 



five components of the spinor as functions of position. In 
these density profiles each vortex (aside from pure phase 
vortices with g = id) will be ornamented by radiating 
density ripples as illustrated in Fig. [51 For example, ac- 
cording to Eq. (|88p . the density of atoms with F^ = m 
is given by 



n(r, 0, m) 



= \Dm2{R)\h;f{r)e^'' + D,n,-i{R)\h9i.r)\ 



3 

b„, cos 



.)■ 



(93) 



where a™ , bm are constants outside the vortex cores, since 
/(r) and g{r) approach 1. While a vortex in a conden- 
sate of a single type of atom does not show any den- 
sity modulation (unless the condensate interferes with a 
second condensate, see e.g. Ref. [501), angular density 
ripples do result for a spinor vortex as a result of the in- 
terference between the / and g components of the spinor 
produced by of the unitary transformation changing the 
quantization axis from the vortex's rotation axis h' to 
the magnetic field direction. If the n' axis happens to 
line up exactly with the axis of the Stern-Gerlach field, 
then there are no radial "interference fringes," but only 
the axes of point vortices with aligned charges will tend 
to line up with B. This is illustrated by the qualitative 
wave function in Section [1111 (See Eqs. ([CTl) . ([HI).) 

Both the order three and order two vortices will be 
visible based on the images of the five spin components. 
One can determine the types of the vortices and their 
axes n' (which are encoded in D{R)) from the average 
magnitude of the densities am together with the ampli- 
tudes h„i and offsets (pm of the density modulations. (An 
order 2 vortex will have cos 2(j) and sin 2(j) Fourier modes 
in addition to the terms given in Eq. ([93l) .) A possi- 
ble difficulty with this method arises because, once the 
five spin components are separated in space, the density 
oscillations in each of them are no longer stable. The 
ensuing dynamics in the clouds could mix the atoms up. 

Distinguishing among vortices with the same rotation 
but different phase winding numbers 9 is not possible 
with this method without resolving the cores. For exam- 
ple, the vortices (i?, ^) and (i?, — %-) have the same den- 
sity patterns, since they differ only by an overall phase 

One would also hope to check some predictions about 
the size and charges of the bound states. One can se- 
lect clusters of vortices in an image of the condensate (if 
there are not too many vortices) and use the methods 
just discussed to identify the vortex charges and check 
that each cluster satisfies Eq. (|54p . Additionally, a sign 
that the vortex clusters are actually bound states is that 
the bound state size depends in the right way on the mag- 
netic field. Now atoms whose ground state is cyclic may 
be difficult to find (®^Rb is fikely to be polar[3^, though 
it may be possible to adjust the interaction parameters 
by applying light fields.). The general considerations of 
this article also apply to spin 3 condensates (see Ref. 
[5l|), as well as to spin 1 condensates and pseudospin -^ 
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condensates as already studied by [ll|, [iJl . So we will es- 
timate the bound state sizes for the more generic phases 
as well as for the cyclic phase. The cyclic phase is unique 
because the misalignment energy is a second order effect 
(see Sec. Ill Ap . For other phases where the misalignment 
energy is a first order effect of the Zceman energy, the 
effective potential would be on the order of quQ. Hence 



La 



for phases with Veff oc q. 



(94) 



In contrast, for the cyclic phase, Eq. ((66)) can be ex- 
pressed in terms of the scattering lengths as 



Lq ^ y/noAa for cyclic phase 

mq 



(95) 



the size of the molecules in the cyclic 
and the size of molecules in 



Since q = ff — 

phase is proportional to 

phases with Veff oc g is proportional to -g. 

The size of the condensate must be large enough hold 
an entire vortex molecule. Substituting BMax from Eq. 
([87)1 into Eqs. dM]) and ((95]) one finds that Lq{BMax) ~ 
—7^ — ~ Cc (for either phase). This size is the magnetic 
healing length of the condensates (and the size of a vortex 
core) and is on the order of 1 /^m. The condensate should 
be narrow in one direction (so that the behavior of the 
order parameter is two-dimensional) but at least several 
times wider than the magnetic healing length in the other 
two directions; in order to measure the field dependence 
of the molecule sizes, one should be able to decrease the 
magnetic field by some factor below BMax without the 
molecules leaving the condensate. 

As a side-comment, vortex molecules probably undergo 
transitions at fields close to B^iax (see Fig. [TO)) since the 
component vortices overlap at larger fields. Absolutely 
stable vortex molecules, like example 3 in Table HI will 
be compressed so that the cores coincide and the vortex 
becomes rotationally symmetric at a finite field. Once the 
components' cores overlap a little bit, being slightly offset 
might not lead to significant savings in kinetic energy. 
On the other hand, when the vortices in a metastable 
molecule are squeezed together, they form an unstable 
tetrahedral vortex. Metastability occurs only when the 
"Coulomb" force keeps the component vortices apart. 
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FIG. 10: The evolution of vortex molecules as the magnetic 
field is increased. The three curves illustrate how the sizes of 
the molecules from Table U for c > 4, might change as the 
strength of the magnetic field is increased. The sizes decrease 
as -J^ . At a certain magnetic field, an absolutely stable vortex 
will become rotationally symmetric and the component cores 
will coincide (molecule 3). Metastable vortices will become 
unstable when a certain magnetic field is reached, indicated 
by the x's terminating the curves corresponding to molecules 
1 and 4. 



More accurate calculations of the vortex fields and en- 
ergies could address other interesting questions. As is 
well-known, unlike the Coulomb interaction between vor- 
tices in an ordinary scalar superfluid, the interaction en- 
ergy cannot be written as a sum of two- vortex interaction 
terms, as indicated by our estimate in terms of the energy 
index. A more accurate understanding of the kinetic en- 
ergy landscape might show that the vortices in a molecule 
can have several spatial arrangements in the core. An- 
other problem that requires more detailed calculations of 
vortex fields is determining whether a vortex with charge 
{R, — ^) is stable: the energy index estimate shows that 
it can break up into vortices with only a finite change in 
energy, but whether the energy increases or decreases is 
not clear yet. 
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VI. CONCLUSION 

We have shown that vortex molecules can be under- 
stood reasonably well based only on simple group the- 
ory and rough energy estimates. Some of these vor- 
tex molecules are actually metastable, and we can study 
their possible break-up "channels," reminiscent of some 
of the decay processes in nuclear physics. (In practice, 
the molecules will probably escape through the surface 
of the condensate before any kind of "ultracold fusion" 
can happen!) 



APPENDIX A: APPENDIX: FINDING THE 
EFFECTIVE ACTION 



Eq. ()39|) is not as difficult to minimize as it appears, 
because of the special symmetry of X2- We must substi- 
tute ip — \/n^X'2. + ^'^1 where 



H - dx2 + aF.^X2 + hFyX2 + cF^X2 + (e + ij)x2t (Al) 

into the energy, Eq. ([55)1 . The perturbation 5i\) is the 
deformation of the tetrahedron, measured relative to its 
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body axes. Let us figure out how many powers of the 
coefficients a,b, . . . to keep at each stage of calculating 
V. It helps to complete the square in Eq. p9|) to get 



VtotW = la{4,^ - mf + i/3(Vi^F^)2 + ^j\4>U\' 



1 



q y^ cos ai cos aj 



4j^F,F,4>, 



(A2) 



where we note that the chemical potential for the cyclic 
state is /i = ano and define 7 — c(3. We also use -^1,2,3 
to stand for F^^y^z- We need to find the minimum of this 
energy only to quadratic order in q. At the end we will 
find that a,b,c, . . . are each linear in q. Since each of the 
squared quantities in V vanishes when a, 6, c, • • • = 0, just 
the linear contributions from a,b,c, . . . give the potential 
to quadratic order in q. The quadratic Zeeman term, 
since it is multiplied by q, also is not needed beyond 
linear order in a, b, c, d, e, /. 

Next find which matrix elements of X2 need to be cal- 
culated to evaluate all these contributions to the en- 
ergy. Only the cross terms between the unperturbed 
part ^/noX2 and the perturbation give linear functions 
of a, b, c, d, e, /. For example, one cross-term contained 
in the quadratic Zeeman contribution is 



i/j'^F^Fyi/j w noxlF.j;FyX2 
+ 2y/7^dixlF,Fy[{aF, 



bFy + cFz + d)x2 



+ {e + iI)X2t]. 



(A3) 



Expanding this gives a sum of matrix elements such 
as X2PxFyX2 and X2tPxFyFzX2- We need only the 
matrix elements of products of at most three i^'s. 
Many of these (e.g., X2FiFjFuX2 when i,j,k are not 
all different, and Xj2tFiFjX2 when i and j are differ- 
ent) are equal to zero because of the 180° symme- 
tries of X2 around the coordinate axes. The numer- 
ical values of the few remaining ones can be worked 
out quickly. Using these matrix elements to calcu- 
late all the terms in Eq. (|A2p produces an expression 
Vtot(a, 5, c, d, e, /, cosai,cosQ;2,cosQ;3). Along the way, 
one notices that each of the variables a, 6, c, . . . con- 
tributes to only one term in the g = potential (the 
first line of Eq. (|A2[) ). The variables a^b,c determine 
the magnetization, d determines the density perturbation 
and e and / determine the singlet-amplitude 0. E.g., 



n = riQ + 2^/n^d 
A4 = 4^/n^a 
^[e\ = 2./T^e. 



(A4) 



Finally, minimize the potential. It can be written as 
a sum of independent quadratic functions of a,6, c, d, e, 



and /: 



V — (— -aJT-o + 2i3'no) + 2anod'^ 



Ay/n^qd 



'8f3no{a^ 



^) + 4\/3noq{a cos a2 cos a^ 



+ b cos ai cos 03 + c cos ai cos 012 ) 
- 27710 (e2 + /2) 

-I- 2q^yno[e{cos^ ai -|- cos^ a2 — 2 cos^ as) 
+ \/3/(cos^ai -cos^a2)]. (A5) 



(Note that the second term, 2qnQ, is the first-order con- 
tribution of the Zeeman energy, which is independent of 
orientation.) 



each quadratic (which gives a = 
- ^y , cosa2 cosaa, . . . ) and combining the results to- 



Minimizing 



gether with the help of Eq. ([57]) gives 



v; 



e// 



2qno - 2 \ — - 



7 



4/3 



+ ( 



3g^_ V 

4/3 7 



)(cosq;^ 



cos a. 



cosa;^); 



(A6) 



all the constant terms can be dropped to give Eq. ([4]). 
Note that the magnetization varies with the orienta- 
tion of the tetrahedron (as can be checked by substitut- 
ing the optimal values for a, 6, c into the magnetization, 
Eq. (jA4|l ). In particular, the c > 4 ground state with 
cosai,cosQ;2,cosQ;3 = ±-7^ has a small magnetization, 
m = =p— ^(1, 1,1); since this has been calculated from 

ip, it is the magnetization relative to the body axes tetra- 
hedron. Comparing this to the magnetic field relative 
to the body axes, i?(cosai,cosa2, cosaa) — -%{^AA): 
shows that the state will become magnetized either par- 
allel or antiparallel to the magnetic field. 

The point of the effective potential is that it al- 
lows us to eliminate the 6 most rigid degrees of free- 
dom corresponding to a,6, c, d, e, and /; then it is eas- 
ier to understand vortex textures by concentrating on 
the rephasing and rotation angles as a function of po- 
sition. The wavefunction in Eq. pip can be parame- 
terized in terms of Euler angles for the rotation, e.g., 
"0 = y/n^e^'^^'' e^^^" e^P^" e^^ X2- The first angle, a, does 
not come into the angles a^ that describe the orienta- 
tion of the magnetic field [Bz) relative to the tetrahe- 
dron. (The tetrahedron can be rotated around the z- 
axis without changing these angles.) One can check that 
cosQi = sin T sin /9, cos a2 = — sinrcosp, and cosaa = 
cos r. Now, working out the kinetic energy and combin- 
ing it with the effective potential gives the "phase-and- 
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rotation-only" energy functional 



APPENDIX B: APPENDIX: 
NONCOMMUTATIVITY OF VORTEX CHARGES 



E.jf = J J d'r^[2iVpr + 2{Vrf + 2{Vaf + 

{Vef -l-4cosTVcr • Vp] 
3q2 
+ (c — 4) — - fees'* T + sin"* r (sin"* p + eos'' p)] . 



(A7) 



This can be solved (in principle) to give the textures 
around sets of vortices and the relative positions of the 
vortices in equilibrium. Each vortex type implies a cer- 
tain type of discontinuity in the four angles as the core is 
encircled. This expression does not seem too easy to use, 
but at least it shows just the two effects we have been 
balancing against one another (kinetic and anisotropy 
energy). The size of a vortex molecule can be esti- 
mated by assuming that the two terms are comparable, 

nph 
mL'i 



c — 4||^. If distances are rescaled by Lq, we then 

find that the energy function has a form that depends 
only on the sign of c — 4. Therefore, in a molecule with 
three vortices, the angles of the triangle they form will 
be independent of all the parameters, including c, even 
though it is dimensionless. 

Eq. (jA7P is derived from Eq. (|32p by determining 
how the tetrahedra are distorted by the quadratic Zee- 
man effect. But the kinetic energy also causes distortions 
of the wave function from the perfect tetrahedral forms, 
and it seems possible that these distortions could lead 
to kinetic effects in the anisotropy term and anisotropy 
effects in the kinetic energy term. However, at the lowest 
order, treating the two terms independently seems cor- 
rect. A simple argument for this (neglecting the kinetic 
energy when finding the anisotropy potential) is that the 
distortion due to the Zeeman term is linear in q while 
the distortion due to the kinetic energy is quadratic in 
q. To see this, think of an ordinary scalar vortex, where 

the density varies as no(l — %) far from the coreQ- The 

2 

amount of "distortion" is — fno- In the cyclic state, 
distortion (i.e., perturbations to the spinor components 
that take it out oi Ai) implies changes in the magneti- 
zation as well as the density. But we may assume that 

2 

these distortions are still of order -j-no- The majority of 
the "pulp" in a molecule's core consists of points whose 
distance is of order Lq from the actual vortex cores (the 
"seeds"), so the amount of distortion can be found by 
substituting r = Lq from Eq. PH)) . Using the relation 
between ac and f3 shows that the fractional distortion 

2 

in the "pulp" regions is of order ^uq ~ (;r0)^; to be 
compared with the deformations of order ^^ that result 
from minimizing Eq. 



To give a complete description of charge conservation 
when the charges are described by the noncommuting 
rotations of a tetrahedron, one needs to give a rule for 
how to multiply the charges of a set of vortices together 
to get the net charge. A convention we used is to multiply 
the topological charges together in order of increasing x- 
coordinates. 

It seems that this definition has an awkward conse- 
quence: does the net vortex charge jump suddenly when 
two of the vortices are reordered, because of the non- 
commutativity of the group of charges? In Fig. [11] vor- 
tices 1 and 2 are interchanged between frames a) and c), 
which suggests that the net charge changes from r3r2ri 
to r3rir2. But this deduction is incorrect, and the net 
vortex charge is actually conserved. 

The resolution of the paradox has to do with the fact 
that the charge of a vortex can only be determined up 
to conjugacy, unless one introduces a systematic conven- 
tion. For example, the charge of a 120°-rotation vortex 
is ambiguous-the rotation could be either P, Q, R, or S, 
and there is no way to distinguish between these because 
the four vertices of the tetrahedra are indistinguishable. 
(Abstractly speaking, the four rotations are conjugate el- 
ements of the group.) In order to identify the charge of 
each vortex, we have to choose a routine for labelling the 
vertices of the tetrahedra nearby. Here is a convention 
that is consistent with the rule for ordering the vortex 
charges. Take a point O far below all the vortices in the 
system and connect it with lines to points just below the 
vortices (see Fig. [TTk). Now identify the base tetrahe- 
dron at O with the standard tetrahedron in Fig. [5^, mak- 
ing a choice from among the twelve possible ways. The 
labelling at O can be communicated to the tetrahedron 
at the end-point of each line, by copying the labelling 
from O to a nearby tetrahedron on the line, and then 
continuing to copy the labelling until the end of the line 
is reached. Now the charges of the vortices can be iden- 
tified by using the labelling of the nearby tetrahedron to 
assign a letter to the rotation axis. 

Now that we have a consistent convention for assigning 
vortex charges, we can show that the net charge of a set 
of vortices does not change when two of them are inter- 
changed. The trick is that the charges of the individual 
vortices do change in such a way that the product charge 
does not change! Between Fig. [TTk and Fig. [TTk . vortex 
2 is moved over vortex 1. Because vortex 2's tether gets 
tangled up with vortex 1 when vortex 1 passes below it, 
its charge gets redefined, as r^"'^r2ri. The other two vor- 
tices' charges do not change. The net charge, obtained 
by multiplying the vortex charges from left to right, is 

r3ri(rr'r2ri) = r3r2ri. 

Thus the net charge does not change. On the other hand, 
there is a sudden jump in the charge of vortex 2, but this 
does not mean that the fields of tetrahedra are changing 
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FIG. 11: a) The convention for assigning vortex charges. For a set of vortices, tethers are drawn directly from the origin to 
an anchor just below each vortex. As long as the tethers are moved continuously, the correspondence between the tetrahedral 
state at the anchor and the standard tetrahedral state does not change. b,c) These show what happens when vortex 1 is moved 
below vortex 2. The dashed lines in a) and the dash-dot-dash lines in c) are the tethers before and after vortex 1 moves. Part b) 
focuses on the tether of vortex 2 , showing how the original tether gets pushed to the side by vortex 1 and is replaced by a new 
tether. Continuing the labelling of the tetrahedron vertices around the original path changes the labelling of the tetrahedron 
just below vortex 2. Hence the charge of vortex 2 is identified differently in c). 



suddenly; the vortex has just been reclassified, with a 
vortex charge that is conjugate to the original charge. 

Here is an interesting consequence of the noncomniut- 
ing charges: the force between a pair of vortices changes 
from repulsive to attractive if a third vortex wanders be- 
tween them. In Fig. [T^ one vortex (the P one) cat- 
alyzes a reaction without touching the other two vor- 
tices involved. (Assume the phases are for the two 
B vortices and ^ for the P vortex.) To figure out 
what happens, keep track of when a vortex's connec- 
tion to the reference point (below the figure) is inter- 
rupted by another vortex. The charge of the vortex 
passing underneath is not changed, and the charge of 
the vortex on top changes to keep the total charge the 
same. This information is sufficient for working out all 
the charges: When P passes below the B on the right, 
the latter vortex changes to a B' = P~^BP, so that the 
net charge is still the same, even though P has moved. 
(This can be used to work out the charges: the net 
charge of the two vortices which have switched has to 
be the same, so PB' = BP so B' = P'^BP.) Next, 
when P passes above the B on the left, the former vor- 
tex changes to a BPB^^ vortex. Now the SU2 matri- 
ces for B and P are —iuy and ~»°"^ "^v '"'^ (using the 
axes associated with X2)- Multiplying out the charges 
shows that B' = A^^ . The force between the original 
pair of vortices (say the P is far away at the begin- 
ning and end) is 1^{Ie{B\Q) - Ie{B,Q) - Ie{BM 
and the force between the vortices they turn into is 



ii^{lE{BA-^),{)) - Ie{A-\{)) - Ie{B,0)). [L is the 
distance between the vortices.) The vortices repel each 
other at first, but after P's intervention, they attract each 
other, as one sees by checking that BA^^ = C, and that 
B'^ is a 27r rotation with energy index 2 while the other 
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FIG. 12: Catalysis by conjugation; the vortex on the right 
moves between the other two vortices changing their repulsion 
to attraction. 



reacting charges are all tt rotations, with index ^. 

We have been assuming that 9 = 0, but a similar reac- 
tion could also happen in the core of a composite vortex 
when g ^ 0; that is why one has to check all the possible 
ways for vortices to wander between one another before 
coalescing or before dividing into groups. 
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This is a small oversimplification. Decreasing the size of 
the space M also stabilizes vortices that would be unsta- 
ble at g = 0-removing a portion of A4 can leave a new 
hole for a vortex-circuit to wrap around. 
The a's for this vortex may be calculated by noting that, 
relative to the tetrahedron, the z-axis of space rotates 
180° about A so that (cos ai, cos a2, cos 0:3), its coordi- 
nates relative to the body axes of the tetrahedron are 
easy to work out as functions of 0. 

Hence, when c > 4, the fluid is ferromagnetic-i.e., it 
becomes spontaneously magnetized either parallel or an- 
tiparallel to the external magnetic field with an m that 
can be calculated from the expressions in Appendix \K\ 
This magnetization spontaneously breaks a symmetry 
since the quadratic Zeeman effect is invariant under 
Fz -^ —Fz- For zero net magnetization, domains will 
form that have magnetization in both directions. Even 
in the presence of a small cubic Zeeman term, one can 
show that there is a discontinuous phase transition (as 
the thermodynamical variable conjugate to magnetiza- 
tion is varied) although the symmetry between ±£ is not 
exact . 
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[55] 



Using the equations, the quadratic Zeeman term distorts 
the ground state y^noxs with increasing magnetic field 
B into ipsiB) = VH^(X(B),0,0,y(B),0)^ with X > 

^,Y < w I because the quadratic Zeeman field favors 
large values of | F^ | . The ground state y^noXa distorts into 
MB) = y^{X'{B),0, -iY'{B),0, X'{B)) with X' > i 
and y < -^. The deformed state ips{B) has a nonzero 
m and a vanishing 6 and V'zC-B) is the other way around. 
Incidentally, when B becomes large enough, the ground 
state will become more symmetrical, with X and X' = 1 
(compare with Ref. ;2&| l . 



[56] The break-up only lowers the net energy by a finite 



amount, about tt- 



[57] 



In 



■g and the fragments interact 

with a short-range repulsion, ex Inr/r^ perhaps; when 
the energy index of the fragments actually decreases the 
force is oc i. Whether (7?,—^) is stable or not we do 
not know, and the answer may depend on c. 
Hence, the direction of the magnetization, n, helps deter- 
mine the orientation of the tetrahedra near each vortex 
core, in spite of Eq. (|46|l. If there are several vortices 
nearby one can try to guess how the tetrahedron fields 
around them fit together. 



